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SECTION  I 


INTRODUCTION 


The  phenomenon  of  combustion  instability  has  plagued  all  types  of  high  heat  release 
combustors  from  industrial  furnaces  to  rocket  engines.  In  general,  however,  the  problem  has 
been  most  severe  in  (light  propulsion  systems  such  as  turbojets,  ramjets,  and  rockets  where 
weight  considerations  dictate  highly  efficient  structures.  In  such  applications,  the  pressure, 
vibration,  and  heat  loads  resulting  from  combustion  instability,  superimposed  on  the  normal 
loading,  are  usually  destructive. 

In  airbreathing  engines,  high  frequency  combustion  instability  problems  were  first 
encountered  in  the  early  1950's,  and  solutions  were  sought  through  mathematical  modeling  and 
analytical  studies  directed  toward  an  understanding  of  the  phenomenon  Unfortunately,  the 
computer  technology  and  analytical  techniques  of  20  years  ago  proved  inadequate,  and  “cut 
and  try"  empirical  approaches  involving  changes  in  flameholders.  combustion  chamber  shape, 
fuel  injection,  ve'oeitv  profiles,  and  flame  piloting  were  attempted.  Fuel  additives  and 
combustion  chamber  baffles  were  also  tested.  Although  some  of  these  approaches,  notably 
baffles,  produced  marginal  improvement,  the  problem  was  not  solved  until  damping  devices  in 
the  form  of  acoustical  absorbers  (screech  liners)  were  introduced.  Screech  liners  are  used 
routinely,  and  high  frequency  instability  is  no  longer  regarded  as  a  problem. 

Larger  and  more  powerful  turbopropulsion  systems  are  presently  being  designed  and 
developed.  Because  of  the  large  physical  dimensions  of  augmentors  used  in  these  systems,  their 
natural  acoustic  modes  have  correspondingly  long  wavelengths;  therefore,  combustion  in¬ 
stability  can  occur  at  very  low  frequencies,  i.e.,  approximately  200  Hz  or. less. 

The  occurrence  of  instability  at  lower  frequencies  makes  use  of  screech  liners  of 
conventional  design  difficult.  To  obtain  adequate  damping,  the  absorbing  devices  are  designed 
so  that  the  resonant  frequency  corresponds  to  the  frequency  of  the  expected  mode  of 
instability.  The  required  cavity  volume  is  inversely  proportional  to  the  square  of  the  resonant 
frequency;  therefore,  low  frequencies  require  large  volumes.  Large  cavity  volumes  can  be 
accommodated  by  increasing  the  augmentor  envelope,  but  this  produce,  an  unacceptable 
increase  in  engine  weight. 

The  problem  of  these  very  low  frequencies,  called  rumble,  have  been  reduced  through 
combined  experimental  and  analytical  techniques.  This  experience  has  emphasized  the  neces¬ 
sity  to  understand  the  fundamental  mechanisms  involved  in  order  to  formulate  a  meaningful 
analytical  effort  and  the  necessity  to  relate  this  effort  to  physical  hardware  and  processes. 

The  main  feature  which  limits  the  usefulness  of  many  combustion  stability  models  is  the 
absence  of  a  direct  correlation  between  the  physical  hardware  ot  a  real  augmentor  and 
combustion  stability.  Experience  has  shown  that  relatively  subtle  altera  ons  in  flameholder 
geometry  can  produce  profound  changes  in  rumble  limits.  For  a  usable  design  and  evaluation 
tool,  the  tnodel  must  he  able  to  relate  directly  to  such  geometry  changes. 
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Such  a  model  was  developed  by  Pratt  &  Whitney  Aircraft  under  two  previous  Air  Force 
contracts.  These  programs  were: 

(1)  Flameholder  Combustion  Instability  Study,  F33615-76-C-2023 

(2)  Lo-Frequency  Augmentor  Instability  Study,  F33615-76-C-2024. 

The  model  which  resulted  from  these  programs  correctly  predicted  the  major  qualitative 
trends  of  augmentor  rumble.  There  were,  however,  some  areas  of  quantitative  improvement 
desired  and  also  some  extensions  of  program  capability. 

The  goals  of  this  effort  were  to  quantitatively  improve  the  predictions  of  the  model  in  the 
areas  of  combustion  efficiency,  rumble  frequency  and  rumble  limits.  Also,  this  program  added 
to  the  model  the  analytical  capability  to  predict  the  augmentor  fuel-air  ratio  distribution  from 
the  input  injector  geometry  and  flow  conditions. 

The  resultant  analytical  models  have  been  combined  into  one  computer  program  which 
has  three  major  modules  as  follows: 

•  The  fuel  injection  module  predicts  the  transient  or  steady-state  behavior  of 
an  input  fuel  injection  system  of  sprayrings  or  spraybars.  The  major  output 
is  a  fuel  flowrate  distribution  in  an  augmentor. 

•  The  combustion  module  predicts  the  efficiency  of  a  turbofan  augmentor  of 
known  geometry  and  inlet  flow  conditions.  1  he  fuel-air  ratio  distribution 
may  be  input  or  calculated  from  the  output  of  the  fuel  injection  module. 

The  major  outputs  are  the  predicted  efficiency  and  the  dynamic  efficiency 
influence  coefficients. 

•  The  rumble  module  predicts  the  stability  of  the  augmentor  at  known 
geometry  and  flow  conditions.  This  module  may  use  either  specified  com¬ 
bustion  data  or  the  results  of  the  combustion  module. 

This  final  report  is  designed  to  provide  a  comprehensive  description  of  the  analytical 
formulation  of  all  three  modules  in  one  document.  Therefore,  it  contains  information  from  the 
final  reports  of  the  previous  studies  as  well  as  the  results  of  this  program.  This  report  is  a 
companion  to  the  User’s  Manual  for  the  Augmentor  .Stability  Management  Program 
(AFAPI.TR-81  -  ). 
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SECTION  II 


TECHNICAL  DISCUSSION 


1.  PROGRAM  ORGANIZATION 

The  Augmentor  Stability  Management  Program  was  structured  as  a  continuation  of  two 
previous  Air  Force-sponsored  programs  which  defined  the  fundamental  models  for  flame 
stabilization  and  low-frequency  stability  in  turbofan  augmentors.  There  were  two  major  goals 
to  the  current  program.  These  were: 

a.  Correct  known  deficiencies  in  the  existing  models  for  flame  stabilization 
and  system  stability. 

b.  Extend  the  capabilities  of  the  models  by  including  an  analysis  for  the  fuel 
injection  system  which  is  capable  of  predicting  the  augmentor  fuel-air  ratio 
as  a  function  of  time. 

The  program  resulted  in  a  computer  program  which  has  the  capability  of  predicting  the 
performance  and  stability  of  a  given  turbofan  augmentor  configuration  at  a  specified  operating 
condition.  The  computer  program  incorporates  the  three  major  models  of  the  analysis  which 
are: 


a.  Augmentor  Fuel  Injection  Model 

b.  Augmentor  Combustion  Model 

c.  Augmentor  Rumble  Model 

The  fundamental  operational  mode  of  the  total  program  is  as  follows: 

•  Fuel  injection  model  —  for  a  specified  fuel  system  geometry  and  specified 
augmentor  conditions,  the  model  predicts  the  fuel  injection  distribution  as  a 
function  of  time.  A  transfer  subroutine  converts  this  into  the  fuel-air  ratio 
distribution  approaching  a  user  specified  flameholder  configuration.  This  is 
the  primary  output  of  this  model. 

•  Augmentor  combustion  model  —  taking  the  now  known  fuel-air  ratio 
distribution  and  input  flameholder  and  augmentor  geometries,  this  model 
will  evaluate  the  level  of  combustion  efficiency  at  specified  augmentor  inlet 
conditions.  The  model  will  also  predict  if  any  of  the  flameholder  elements 
are  beyond  the  local  stabilization  limit.  The  model  then  generates  a  series  of 
normalized  combustion  gain  factors  which  are  to  be  used  in  the  rumble 
analysis. 

•  Augmentor  rumble  model  —  using  the  given  geometry,  conditions,  efficiency 
and  gain  factors,  this  model  predicts  the  dynamic  gain  and  phase  angle  of 
low-frequency  oscillations.  An  instability  exists  whenever  a  gain  of  unity  or 
more  is  shown  with  a  phase  angle  of  zero  degrees. 

The  overall  computer  code  includes  all  three  of  these  models.  They  may,  however,  be  run 
in  any  separate  or  sequential  fashion  if  the  user  specifies  the  required  transfer  data. 
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B.  TECHNICAL  RESULTS 


1.  Augmentor  Fuel  Injection  Model 

a.  General 

Aircraft  gas  turbine  engines  with  thrust  augmentation  share  a  common  problem  of 
requiring  a  fuel  delivery  system  in  a  high-temperature  environment.  During  nonaugmented 
operation,  augmentor  sprayrings  or  spraybars  soak  in  the  turbine  exhaust  gas  and  are  heated 
to  the  gas  temperature.  During  augmentation  the  fuel  system  is  cooled  by  the  (lowing  fuel. 
Problems  can  occur  during  the  transition  from  nonaugmented  to  augmented  operation.  The 
transient  start  process  requires  the  filling  of  hot  fuei  lines  with  fuel,  and  may  result  in  fuel 
boiling,  surging,  and  failure  to  achieve  and  sustain  augmentor  ignition.  Depending  on  the 
particular  fuel  system  configuration,  different  engines  may  have  significantly  different  prob¬ 
lems  related  to  the  transient  start  of  augmentation. 

b.  Physical  System 

An  augmentor  fuel  system,  as  defined  in  this  model,  consists  of  a  fuel  control  system,  a 
fuel  distribution  system,  and  necessary  connecting  plumbing.  Two  types  of  distribution 
systems  are  considered  in  the  model,  circumferential  sprayrings  and  radial  spraybars,  and  both 
types  will  henceforth  be  referred  to  generically  as  “injectors.”  A  typical  sprayring,  shown  in 
Figure  1,  is  supplied  with  fuel  at  one  or  more  inlet  locations  along  its  .  ircumference.  Fuel  spray 
jets  may  be  either  constant  or  variable  area  orifices,  and  are  distributed  ah>ng  the  sprayring.  In 
a  sprayring  fuel  system,  multiple  rings  of  different  diameters  are  positioned  coaxially  in  the 
augmentor.  and  may  be  supplied  with  fuel  individually  or  in  various  multiple  combinations. 


In  a  spraybar  system,  radial  spraybars  are  positioned  circumferentially  around  the 
augmentor.  Fuel  supply  is  from  one  end,  and  spraybars  may  be  jointly  supplied  in  multiple 
combinations  or  zones. 

Fuel  control  systems  operate  to  control  either  fuel  inlet  How  rate,  system  pressure,  or  a 
combination  of  the  two.  Furthermore,  provisions  may  be  made  for  rapid  system  fill  during 
start  transient,  or  for  sequenced  starting  of  multiple  zone  systems. 

During  nonaugmented  operation  the  augmentor  fuel  supply  system  is  filled  with  air.  and 
soaks  in  the  hot  turbine  exhaust  environment.  During  the  transition  to  augmented  operation, 
fuel  fills  the  system,  displacing  the  air  and  cooling  the  heated  injector.  The  incoming  fuel  is 
heated  and  boils  as  it  contacts  the  hot  injector  walls.  The  air  and  fuel  vapor  inside  the  injector 
are  compressed  by  the  incoming  fuel,  causing  system  pressure  to  rise,  which  in  turn  affects 
both  the  heat  transfer  to  the  fuel,  and  the  rates  of  fuel  flow  and  air  flow  out  of  the  system.  The 
transient  phenomena  associated  with  start  of  augmentation  are  controlled  by  the  complex 
interaction  of  heat  transfer,  single  and  two-phase  flow  of  air  and  fuel  out  of  the  system,  and 
thermodynamic  behavior  of  the  system  components. 


Upper  Inlet 


Figure  1 .  Schematic  of  Sprayring  During  Fill  Process 


Several  simplifications  are  made  in  describing  the  physical  system  by  an  analytical 
model.  The  three  most  signifcant  assumptions  are  the  following: 

1.  The  fluid  system  can  be  described  as  three  regimes:  a:  liquid  fuel,  and 
two-phase  fuel  (liquid-vapor).  These  regimes  do  not  mix. 

2.  The  two-phase  fuel  regime  can  be  described  as  a  homogeneous  mixture  of 
liquid  and  vapor. 

3.  The  pressure  throughout  the  fuel  system  is  uniform  at  any  time  (i.e., 
pressure  losses  within  the  system  are  small  compared  to  temporal  rates  of 
pressure  change). 
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Figure  2  illustrates  the  mouel  assumptions  as  applied  to  a  spraybar.  The  illustrated 
sprayhar  consists  of  two  sections,  a  cold  section  oriented  out  of  the  hot  turbine  exhaust,  and 
a  heated  section  in  the  turbine  exhaust.  The  fuel  spray  jets  are  assumed  to  be  distributed 
uniformly  over  the  entire  length  of  the  spraybar.  The  coordinates  shown  in  Figure  2  define  the 
length  of  the  subcooled  liquid  regime  (X„),  the  length  of  the  two-phase  regime  (X  -  X„),  and 
the  length  of  the  air  regime  (Za). 
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Figure  2.  Schematic  of  Spraybar  Showing  Fluid  Regimes 


To  begin  the  analysis  of  the  augmentor  start  transient,  the  initial  conditions  must  be 
specified.  Initially,  the  injector  is  assumed  to  be  air-filled  and  the  metal  temperatures  are 
assumed  equal  to  the  local  turbine  exhaust  temperatures.  As  fuel  enters  the  system,  toe  air  is 
compressed  and  flows  out  through  the  spray  jets.  Fuel  is  heated,  eventually  boiling,  and  the 
fuel  vapor  adds  to  pressurization  of  the  system.  Air  is  gradually  purged  from  the  system,  and 
the  transient  continues  until  the  injector  is  cooled  to  the  level  determined  by  steady-state  heat 
transfer  rates. 

c.  Continuity  and  System  Dynamics 

Considering  air,  liquid  fuel,  and  fuel  vapor  as  separate  species,  mass  continuity  equations 
can  be  written  for  each  specie.  In  an  injector  of  cross  sectional  area  A,  the  liquid  fuel  mass  (M) 
is: 


Mt  =  ,>A  [X  -  «(X  -  X„)|  (1) 

where  p  is  the  mean  liquid  density,  «  is  the  vapor  void  fraction,  and  (X,  X,,)  are  regime  lengths 
defined  in  Figure  2.  Similarly,  the  masses  of  vapor  and  air  in  the  system  are: 

Mv  =  py A  <v  (X  -  X„)  (2) 

Ma  =  pa A  Za  (3) 

where  pv  and  pa  are  the  densities  of  vapor  and  air.  Specie  continuity  is  described  by: 


dt 

dM. 

dt 

iM. 

dt 


=  mu  -  m|(,  -  mv 

(liquid) 

(4) 

=  niv  -  mv„ 

(vapor) 

(n) 

=  ma 

(air) 

(B) 

where  (m|,,  ni|„)  are  liquid  flow  rates  in  and  out  of  the  system,  m,  is  the  rate  of  vapor 
generation,  and  (mv,„  m„)  are  rates  of  vapor  and  air  flow  out  of  the  system. 


Combining  Equations  (1)  through  (6),  assuming  air  compression  to  be  isentropic,  and 
vapor  compression  to  be  at  constant  quality,  a  set  of  continuity  equations  is  produced  relating 
system  pressure  and  liquid  and  vapor  volume  Solving  the  set  of  equations  simultaneously  for 
the  interdependent  variables  results  in  a  set  of  equations  defining  the  rates  of  change  of 
pressure,  liquid  volume,  and  vapor  volume  in  the  system.  The  resulting  continuity  set.  written 
for  a  spraybar  fed  from  one  end  is: 
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(9) 


The  dependent  variables  are  system  pressure  (P),  fuel-air  interface  location  (X),  two-phase 
interface  (X„),  and  vapor  void  fraction  («).  From  the  nature  of  this  equation  set,  it  is  apparent 
that  heat  transfer  and  flow  calculations  must  be  performed  prior  to,  or  simultaneously  with, 
evaluation  of  the  set,  and  are  in  turn  dependent  upon  the  results. 


d.  Heal  Transfer  Model 

Because  of  the  strong  dependence  of  system  performance  on  fuel  vaporization  rate  (mv), 
as  evidenced  in  the  continuity  equations,  provision  of  an  accurate  description  of  the  fuel  heat 
transfer  characteristics  is  a  prerequisite  to  the  realistic  simulation  of  the  physical  system. 

Two  principal  heat  transfer  processes  occur:  heat  flows  from  the  injector  to  the  fuel 
during  the  transient  cooling  of  the  system,  and  as  the  injector  cools,  heat  flows  from  the  hot 
turbine  exhaust  gases,  establishing  a  steady-state  process.  The  transient  process  is  the  primary 
source  of  heat  flow  during  system  fill  and  cooldown. 

Heat  transfer  from  the  injector  to  the  flowing  fuel  can  occur  in  one  of  several  modes, 
combining  simple  convection  with  various  boiling  modes.  Because  of  the  high  initial 
temperature  of  the  injector,  film  boiling  of  the  fuel  is  possible.  Because  the  heat  source  is  from 
energy  storage  in  the  injector,  the  source  temperature  is  the  controlling  parameter  in  the  heat 
transfer  process.  As  the  system  cools,  heat  transfer  modes  will  vary  from  forced  convection  film 
boiling,  to  forced  nucleate  boiling,  to  simple  convection.  The  local  instantaneous  difference 
between  the  injector  temperature  and  the  fuel  liquid  saturation  temperature  is  considered  to 
be  the  controlling  parameter  for  determination  of  the  heat  transfer  regime.  Figure  3  shows  the 
principal  differences  in  the  several  heat  transfer  modes.  Although  vapor  is  generated  by  boiling 
in  subcooled  fuel,  the  vapor  is  rapidly  condensed  in  the  colder  fluid  away  from  the  wall.  Vapor 
accumulation  only  begins  when  the  bulk  fluid  temperature  equals  the  liquid  saturation 
temperature  (Ts().  During  saturation  boiling,  the  fuel  bulk  temperature  (T h)  increases  from 
the  saturated  liquid  to  the  saturated  vapor  state  at  the  location  of  total  vaporization,  or 
dryout.  The  various  regimes  are  shown  again  in  Figure  4  as  a  function  of  temperature 
difference  between  wall  and  saturated  liquid  state  (ATW).  Three  specified  temperature  dif¬ 
ferences  are  used  to  identify  heat  transfer  regimes:  the  onset  of  nucleate  boiling  (AT,,),  the 
maximum  nucleate  boiling  condition  (AT*),  and  the  minimum  film  boiling  condition  (ATct)- 
In  each  heat  transfer  regime  the  liquid  fuel  may  be  in  a  thermodynamic  state  ranging  from 
subcooled  to  saturated. 
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Figure  ■'!  Heat  T 'ransfer  Processes  in  Sprayring 
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Figure  /.  Fuel  Heat  Transfer  Regimes 

Figure  5  illustrates  the  thermodynamic  process  path  as  a  blended  hydrocarbon  fuel  (such 
as  JP-5)  is  heated  at  constant  pressure.  The  heat  transfer  characteristics  in  each  heat  transfer 
regime  are  defined  by  an  empirical  model.  Figure  '  shows  the  logic  for  selection  of  the 
appropriate  empirical  model  for  description  of  the  local  heat  transfer  to  the  fuel. 
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Figure  5.  Path  of  Thermodynamic  Process  for  Fuel  Flow  in  Sprayring 
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Figure  6'.  Logic  for  Heat  Transfer  Model 

The  specific  models  used  are  as  follows: 

1.  Nonboiling  convection  is  described  by  the  Dittus-Boelter  relation  for 
turbulent  flow,  with  a  constant  lower  limit  of  the  Nusselt  number  for 
laminar  flow. 

2.  Nucleate  boiling  is  described  by  interpolation  between  the  convection 
model  and  the  maximum  nucleate  boiling  heat  flux,  described  by 
Kaminski'  for  JP-5  by: 

qm«*  =  27.4  +  0.351  ATsr  V  (Btu/hr-ft2)  (10) 

^here 

yi\r  =  Tsl  -  Th,  °R 

V  =  fuel  velocity,  ft/sec 


!i.  Saturated  film  boiling  is  described  by  (ilickstein  and  Whitesides'  lor  light 
hydrocarbons  by: 

u 

Nu  -  0. 10f>8  Re,"-'*-"  Pr,"4  ^  ~ ^  ill, 

where  the  Reynolds  number  ( He,  I  is  based  on  the  local  fuel  velocity  and  p(l 
is  the  bulk  density. 

4.  Subcooled  film  boiling  is  described  for  .IP-5  by  a  modified  Dittus-Hoelter 
equation  of  the  form: 

Nu  ■=  0.529  Re,"8  Pr"4  Sn"8  (12) 

where  Sn  is  the  Sternum  Number,  defined  by 


(14) 


The  vapor  density  (/>,)  is  evaluated  at  the  arithmetic  average  film 
temperature,  and  H,-„  is  the  constant  pressure  latent  heat  of  vaporization. 

5.  Roiling  in  the  transition  regime  is  treated  by  interpolation  between  max¬ 
imum  nucleate  boiling  and  minimum  film  boiling. 

e.  Flow  Model 

Accurate  evaluation  of  the  fuel  and  air  How  rates  out  of  the  injector  is  an  implicit 
requirement  for  the  evaluation  of  the  continuity  equations.  In  addition,  fuel  velocity  within  the 
system  is  required  for  evaluation  of  the  heat  transfer  to  the  fuel.  Since  the  local  internal 
velocity  is  dependent  on  the  distribution  of  overboard  flow  (out  of  the  system),  and  the  local 
thermodynamic  state,  it  is  apparent  that  the  overboard  mass  flow  rates  and  internal  heat 
transfer  are  implicitly  related. 

Characterization  of  the  overboard  flow  is  simplified  by  defining  several  flow  regimes. 
Flow  of  air  from  the  injector  is  described  as  one-dimensional  compressible  flow,  and  may  be 
subsonic  or  choked  flow.  Fuel  flow  regimes  are  somewhat  more  complex,  and  may  be 
distinguished  by  examining  pressure  potential  throughout  the  injector.  If  Psh  is  the  bulk 
saturation  pressure  corresponding  to  the  local  fluid  bulk  temperature,  then  Psb  will  vary 
throughout  the  system.  Figure  7  shows  the  variation  of  Ps),  over  the  length  of  a  heated  injector, 
compared  with  the  system  pressure  (P)  and  the  augmentor  pressure  (P„).  The  regimes  are 
defined,  (I)  Ps|,  <  P,„  (II)  subcooled  with  Pst)  >  P,„  and  (III)  two-phase  flow: 

1.  Regime  I  —  Fuel  is  subcooled  and  vapor  pressure  is  less  than  augmentor 
pressure.  Fuel  flow  is  described  as  liquid  flow  through  an  orifice  with  the 
pressure  difference  (P  P„). 

2.  Regime  II  Fuel  is  subcooled,  but  vapor  pressure  exceeds  augmentor 
pressure.  Fuel  is  assumed  to  vaporize  as  it  flows  through  the  spray  jets, 
limiting  discharge  pressure  to  P.,,.  Flow  is  described  as  orifice  flow  with  the 
driving  pressure  difference  (P  -  Ps|,). 
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3.  Regime  III  -  Two-phase  flow  with  bulk  boiling.  Fuel  flow  is  choked  in  the 
spray  jets,  the  flow  process  being  described  as  isentropic  expansion  from 
system  pressure  (F),  to  a  critical  pressure  (P*).  The  ratio  of  critical  to 
system  pressure  (P*/P)  is  a  function  of  local  vapor  void  fraction,  described 
by  Starkman,'1  et  al.  and  ranges  from  0.72  to  approximately  0.5.  The  vapor 
fraction  of  the  total  flow  in  the  two-phase  regime  is  assumed  to  equal  the 
initial  local  vapor  quality. 
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Figure  7.  Definition  of  Fuel  Flou  Regimes 


f.  Program  Logic 

The  described  models  for  continuity,  heat  transfer,  and  mass  flow  have  been  programmed 
for  computer  evaluation.  Required  thermodynamic  and  thermophysical  properties  of  several 
fuels  have  been  incorporated  in  the  program  as  empirical  correlations.  Although  simultaneous 
evaluation  of  the  various  models  would  be  desirable,  sequential  evaluation  is  used  in  the 
program  because  of  its  inherent  simplicity. 

Transient  simulation  of  the  augmentor  start  sequence  begins  with  a  defined  initial  system 
condition.  Evaluation  of  the  various  models  proceeds  in  a  time-marching  mode,  with  an  update 
of  all  models  and  properties  at  each  time  step.  Spatial  evaluation  of  heat  transfer,  mass  flow 
and  injector  and  fuel  temperature  distributions  are  conducted  with  finite  element  analysis  at 
each  time  step.  The  operational  sequence  of  the  computer  evaluation  is  outlined  in  the  logic 
diagram  shown  in  Figure  8. 

Time  marching  evaluation  continues  until  the  system  becomes  totally  liquid-filled. 
Referring  to  Equation  (7)  for  rate  of  pressure  change,  it  is  noted  that  the  denominator  tends 
to  zero  as  the  system  becomes  totally  free  of  air  and  fuel  vapor.  For  that  case,  the  numerator 
of  Equation  (7)  must  be  identically  zero,  defining  the  incompressible  flow  condition. 
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Figure  S  Flint  Logic  for  Augmentor  Fuel  System  Model 

Continuity  Equation  (7)  defines  rate  of  change  of  system  pressure  as  a  function  of  several 
independent  variables,  including  inlet  fuel  flow.  As  a  result,  the  program  can  most  conve¬ 
niently  be  used  to  simulate  fuel  systems  with  inlet  (low  rate  controlled  as  a  function  of  time. 
The  model  can  also  be  readily  used  for  multiple  zone  systems.  To  allow  evaluation  of  such 
systems,  the  program  has  been  configured  to  operate  in  a  multiplex  mode,  sequentially 
evaluating  at  each  time  step  the  characteristics  of  each  zone.  Operated  in  this  mode,  the 
program  can  simulate  multiple  augmentor  fuel  systems  with  different  zones  programmed  to 
start  in  sequence,  at  the  same  time,  or  in  any  other  combination. 
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g.  Analytical  Results 


The  pregram  has  been  compared  with  engine  test  data  at  several  different  operating 
conditions,  with  good  results.  A  single  sprayring,  with  controlled  inlet  fuel  flow,  was  exercised 
through  the  start  sequence  with  the  engine  operating  both  at  sea-level  static  and  also  at  high 
altitude  subsonic  conditions.  Figure  9  shows  the  inlet  How  rate  schedule  used  in  the  SLS  test, 
and  also  used  as  program  input.  Figure  10  shows  the  resulting  fuel  system  pressures,  measured 
in  the  test  and  predicted  by  the  simulation.  Figure  11  and  12  show  the  same  comparison  for  a 
test  at  Mach  No.  0.8  and  50,000  ft  altitude.  In  both  comparisons  the  model  results  agree  well 
with  the  engine  data,  although  some  differences  are  apparent  in  predicted  pressure  levels. 
Continuing  model  development  is  currently  in  progress  to  refine  the  quantitative  accuracy  of 
the  model. 
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Future  9.  Inlet  Fuel  Flow  in  Sprayring  Test  at  SLS  Conditions 
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Figure  10.  Comparison  of  Measured  and  Predicted  Pressures  in  Fuel  Sys¬ 
tem  During  Augmentor  Start  at  SLS  Condition 


Figure  11  Inlet  Fuel  Flow  in  Sprayring  Test  at  Mach  0.8/50,000  ft  Altitude 
Conditions 
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2.  Development  of  the  Combustion  Model 


a.  Rumble  Mechanism 

The  nature  of  the  feedback  mechanism  which  drives  rumble  is  the  classical  pressure  and 
velocity  response  of  the  heat  release  mechanism  causing  variations  in  the  average  augmentor 
exit  temperature.  These  variations  result  in  oscillations  in  average  nozzle  inlet  pressure 
through  the  choked  flow  relationships  for  a  constant  mass  flowrate.  This  results  in  further 
oscillations  in  the  flow  conditions  around  the  sprayring  and  flameholder  region.  The  relative 
amplitudes  and  phase  angles  are  functions  of  the  geometry  and  operating  conditions  of  the 
augmentor  and  fan  duct. 

'The  analysis  and  model  for  rumble  were  structured  to  rely  as  heavily  as  possible  on 
analyses  of  the  physical  processes  of  the  fuel  preparation.  Ilame  stabilization  and  flame 
propagation.  In  this  manner  it  was  felt  that  maximum  utility  of  the  model  would  be  realized. 
All  too  frequently  combustion  stability  models  are  generated  with  built-in  correlation  cons¬ 
tants  which  are  generic  to  one  form  of  combustor  only.  Whenever  a  geometry  variation  is 
performed  there  is  no  guideline  for  the  required  change  in  those  constants,  and  the  utility  of 
the  model  is  limited  until  sufficient  experience  is  obtained.  In  this  model  format  we  have 
attempted  to  remove  that  restriction. 

The  mechanism  of  response  of  the  heat  release  process  to  variations  in  pressure  and 
velocity  over  the  sprayring  and  flameholder  region  is  described  below  with  reference  to 
Figure  Ilk 

•  A  variation  in  pressure  and  velocity  over  the  fuel  injection  spraybar  results 
in  a  variation  in  the  average  air  mass  flowrate  and  thus  local  fuel-air  ratio, 
since  the  fuel  flowrate  is  essentially  constant. 

•  At  the  flameholder,  the  pressure  and  velocity  oscillation  results  in  a  vari¬ 
ation  in  the  rate  of  formation  of  the  flameholder  wake  fuel-air  ratio  from  the 
available  liquid  phase  fuel. 

•  The  pressure/velocity  oscillation  also  results  in  a  variation  in  the  level  of  the 
wake  reaction  efficiency  at  the  above  level  of  fuel-air  ratio.  This  results  in  a 
variation  in  the  rate  of  initial  flame  propagation  into  the  approach  stream 
fuel-air  mixture. 

•  The  pressure/velocitv  oscillation  also  results  in  a  variation  in  the  rate  of 
transverse  flame  spreading  at  the  local  flame  speed  due  to  the  influences  of 
local  pressure  and  turbulence  on  flame  speed. 

•  The  final  result  is  a  variation  in  the  average  value  of  the  nozzle  inlet 
temperature  due  to  the  combined  effects  of  the  previous  responses. 

This  variation  occurs  at  a  time  delay  equal  to  the  transport  delay  between  the  spravrings  and 
the  nozzle. 

This  process  includes  two  major  mechanisms  which  cause  a  response  in  nozzle  inlet 
temperature  to  oscillations  in  augmentor  inlet  pressure  or  velocity.  The  variation  in  local 
fuel-air  ratio  (i.e.,  axial  variation  at  fixed  overall  level)  causes  an  axial  variation  in  heat  release 
which  is  felt  as  a  temporal  oscillation  in  nozzle  inlet  temperature.  Imposed  onto  this  variation 
is  the  variation  in  axially  local  heat  release,  at  the  value  of  fuel-air  ratio,  caused  by  response 
of  the  flame-holding  process  to  the  pressure  or  velocity  oscillations. 
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ipure  Id.  Two-Phase  Fuel  Ducted  Flame  Model 


b.  Combustion  Model 


The  basic  framework  of  Ihe  model  will  he  described  with  references  to  Figure  hi.  'This 
illustration  may  be  visualized  as  representing  a  two-dimensional  sector  of  a  flameholder  array. 
A  full  augmentin'  analysis  would  require  a  multistreamtube  solution  to  represent  the  full 
flameholder. 

The  augmentor  inlet  conditions  are  known  in  terms  of  area,  pressure,  velocity, 
temperature,  turbulence  intensity  and  vitiation  level  (if  any).  The  required  overall  fuel-air 
ratio  and  thus  fuel  flowrate  from  the  spravbar  is  specified  in  terms  of  flowrate  versus  pressure 
drop  and  droplet  size  and  distribution  also  versus  pressure  drop.  Once  specified,  the  Howrate 
of  fuel  or  overall  fuel-air  ratio  and  duct  static  pressure  will  define  a  droplet  size  distribution. 

At  this  point  of  injection  a  simple  enthalpy  balance  is  performed  to  evaluate  the 
percentage  of  the  fuel  which  vaporizes  due  to  the  adiabatic  throttling  process  of  injection  from 
the  high  pressure  sprayring  into  the  low-pressure  augmentor.  For  this,  the  fuel  properties  and 
fuel  temperature  in  the  sprayring  must  he  known.  The  fuel  flowrate  which  remains  in  the 
liquid  phase  is  placed  into  5  or  10  fuel  droplet  size  groups  which  represent  equal  mass  flowrate 
distributions  of  the  spray  distribution  curve  for  the  particular  sprayring  injector. 


These  fuel  droplets  are  allowed  to  accelerate  with  the  airstream  towards  the  flan  'holder 
with  concurrent  droplet  vaporization.  The  vaporization  analysis  utilizes  a  forced  convection 
model  for  combined  heat  and  mass  transfer.  A  finite  difference  scheme  integrates  the 
accelerating  vaporizing  droplet  lifetime  until  either  the  flameholder  plane  is  reached  or  the 
droplet  diameter  goes  to  zero.  This  is  performed  for  each  droplet  size  group.  The  droplets  are 
treated  as  spherical  for  the  acceleration  and  vaporization  analysis. 

It  is  at  this  location,  the  flameholder  leading  edge,  that  the  combustion  model  deviates 
from  classical  flameholder  or  combustion  analyses.  The  classical  explanation  for  the  mecha¬ 
nism  by  which  a  flameholder  functions  is  that  the  wake  of  the  bluff  body  serves  as  a  volume 
for  the  reaction  of  the  vapor  phase  mainstream  fuel-air  mixture.  Ignition  of  the  high-speed 
mainstream  mixture  is  achieved  by  the  hot  wake  reaction  products  in  the  shear  layers  aft  of 
the  flameholder.  The  degree  of  reaction  achieved  in  fhe  wake  volume  is  dependent  on  the  ratio 
of  average  residence  time  to  average  reaction  rate.  'Fhe  reaction  rate  depends  on  the  entering 
vapor  fuel  air  ratio  and  operating  pressure  and  temperature,  i.e..  kinetic  rate  constants. 

The  problem  is  that  in  a  relatively  cool  airstream,  such  as  the  25(1  deg  fan  duct  exit 
temperature,  and  overall  fuel-air  ratios  which  typically  react  well,  the  degree  of  droplet 
vaporization  is  very  low.  So  low,  in  fact,  that  the  vapor  phase  fuel-air  ratio  at  the  flameholder 
is  well  below  the  lean  flammability  limit  for  -IP-tvpe  fuels.  However,  these  conditions  do 
product1  stable  flames  in  actual  burners.  Obviously,  there  is  some  additional  mechanism  by 
which  vapor  phase  fuel  is  generated  and  mixed  into  the  flameholder  w:ake. 


The  mechanism  which  has  been  identified  for  this  process  is  the  formation  of  a  liquid 
film  on  the  surface  of  the  flameholder  by  impingement  of  the  fuel  droplets.  This  surface  film 
is  partially  vaporized  by  the  heat  flux  from  the  hot  wake  through  the  surface  of  the 
flameholder  into  the  film.  This  vaporized  fuel  enters  the  recirculation  zone  through  the  shear 
layer  and  provides  the  bulk  of  the  vapor  fuel  for  the  wake  reaction  process. 

For  the  purposes  of  the  model,  the  process  of  fuel  collection  on  the  surface  of  the 
flameholder  is  analyzed  as  spherical  droplet  trajectories  through  the  flow  field  as  it  moves 
around  the  flameholder.  The  analysis  is  performed  for  each  of  the  initial  drop  size  groups 
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utilizing  the  diameter  which  remains  after  the  droplet  has  experienced  partial  vaporization 
between  the  spraybar  and  the  flameholder.  This  analysis  defines  the  fuel  collection  rate  onto 
the  surface  of  the  flameholder,  when  integrated  over  the  droplet  size  groups. 

Once  the  film  is  established  on  the  surface,  the  rate  of  vaporization  is  analyzed  by 
assuming  a  Nusselt  number  form  of  forced  mass  transfer  driven  by  the  vapor  pressure  of  the 
film.  A  combined  film  heating  and  vaporization  solution  is  performed  as  a  finite  difference 
analysis  along  the  surface.  The  heat  flux  is  evaluated  from  a  wake  film  coefficient  and  an 
assumed  wake  temperature.  Since  the  wake  temperature  is  a  function  of  the  fuel  vaporization 
rate  and  the  wake  reaction  level,  an  iterative  solution  is  required  between  this  initial  guess  of 
wake  temperature  and  the  calculated  value. 

As  the  finite  difference  solution  proceeds,  a  certain  amount  of  liquid  fuel  accumulates 
which  has  experienced  heating  but  no  vaporization.  This  fuel  is  lost  to  the  near  field 
stabilization  process  but  does  enter  into  the  far  field  reaction  during  turbulent  flame  spread¬ 
ing.  The  energy  required  to  raise  this  fuel  from  its  collection  temperature  to  the  final  film 
temperature  represents  a  heat  loss  from  the  wake  reaction. 

For  the  evaluation  of  the  wake  reaction,  the  wake  is  treated  ac  a  well-stirred  chemical 
reactor  operating  with  gaseous  fuel.  To  perform  this  calculation,  the  volume  and  mass  influx 
rate  must  be  known.  These  values  are  produced  from  published  data  on  the  relative  size  and 
recirculation  rates  behind  bluff  body  stabilizers  as  functions  of  geometry  and  operating 
conditions.  (Specific  source  references  will  be  given  as  the  analyses  are  developed  in  a  later 
section  of  this  report.)  These  data  have  been  reduced  to  a  series  of  empirical  functions  which 
relate  wake  volume  and  recirculation  rate  to  flameholder  blockage,  geometry,  approach  flow 
conditions,  and  turbulence  level. 

The  reaction  efficiency  in  this  known  size  wake  reactor  is  analyzed  assuming  that  the 
process  proceeds  as  a  single  step,  second  order  reaction  with  the  inefficiencies  represented  as 
unreacted  CO  for  lean  operation  or  unreacted  fuel  for  rich  operation.  The  entering  mass 
flowrate  is  balanced  against  the  mass  consumption  rate  at  an  unknown  final  level  of  efficiency. 
A  straight-forward  solution  for  the  reaction  efficiency  and  thus  wake  temperature  proceeds. 
Any  external  heat  gain  or  loss  mechanisms  are  used  to  adjust  the  reaction  rate  and  increase  or 
reduce  the  final  reactor  efficiency. 


At  the  known  level  of  entering  air  and  droplet  vaporized  fuel  flowrate  plus  surface 
vaporized  fuel,  the  iteration  is  performed  as  follows: 

Recirculation  rate  calculated 

Wake  temperature  assumed 

Surface  vaporization  calculated 

W'ake  fuel-air  ratio  thus  known 

Wake  reaction  level  at  this  f/a  calculated 

Wake  temperature  thus  known 

Iterate  back  through  surface  vaporization 

At  this  point  the  wake  conditions  are  known,  and  the  fuel-air  ratio  distribution  around 
the  flameholder  is  known.  The  turbulent  flame  spreading  into  the  free-stream  is  initiated  in 
the  shear  layer  by  the  hot  wake  products.  The  degree  of  perfect  initiation  depends  on  the 
excess  thermal  energy  available  from  these  products.  Flame  initiation  is  a  go  or  no-go 
phenomenon  and  statistical  in  nature,  since  a  lack  of  perfect  initiation  physically  results  in 
local  regions  along  the  surface  of  zero  ignition.  As  the  temperature  of  the  wake  products  is 
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reduced,  the  percentage  of  the  area  which  fails  to  ignite  increases.  As  these  regions  increase  in 
number  and  size,  a  greater  portion  of  the  heal  liberated  in  those  areas  which  do  ignite  is 
transferred  into  adjacent  unlit  regions  rather  than  into  transverse  flame  propagation.  Since  the 
model  is  based  on  uniform  flame  initiation,  this  process  is  observed  as  a  slower  rate  of  flame 
spreading.  To  produce  this  result,  the  model  relates  the  flame  speed  to  the  ideal  turbulent 
value  and  the  wake  efficiency  level. 

The  form  of  the  flame  spreading  model  is  a  thin  sheet  flame  front  propagating  into  a 
fuel  air  mixture  at  a  velocity  which  was  accelerated  by  the  flameholder  blockage.  The  specific- 
analysis  for  the  turbulent  flame  speed  follows  Karlovitz,  where  the  turbulent  flame  speed  is 
related  to  the  laminar  flame  speed,  at  the  approach  conditions  and  the  local  fuel-air  ratio,  and 
the  local  value  of  turbulent  velocity.  This  latter  term  is  evaluated  from  the  turbulence 
generated  by  the  flameholder  and  axially  decayed  in  a  10  L/I)  length,  based  on  the  width  of 
the  effective  jet  flow  between  the  flamehoiders,  as  a  function  proportional  to  x-1  >.  The  final 
value  after  decay  is  the  free-stream  intensity. 

An  additional  term  is  added  to  the  calculated  value  of  the  turbulent  flame  speed  to 
account  for  the  sensitivity  of  the  reaction  rate  to  integrated  efficiency  and  flame 
self-turbulence.  These  terms  are  modeled  as  a  multiplier  whose  value  depends  on  the  local 
efficiency  and  has  a  value  of  I  at  zero  and  100%  efficiencies.  The  peak  value  is  2  at  50% 
efficiency. 

The  transverse  position  of  the  (lame  front  is  found  by  a  finite  difference  integration  of 
the  local  flame  front  into  the  approach  flow.  Due  to  the  sustaining  effect  of  the  local  heat 
release,  the  approach  velocity  is  axially  retained  at  the  value  accelerated  by  the  flameholder. 


The  overall  efficiency  is  found  at  the  nozzle  entrance  through  integration  of  the  flame 
spreading  heat  release  evaluated  from  the  preceding  analyses. 

1'he  combustion  analyses  are  computerized  into  one  cohesive  pro-  am  which  combines 
the  models  for  the  various  processes.  The  computerization  approach  selected  is  one  which 
utilizes  separate  subroutines  for  each  of  the  process  sections  of  the  model.  Although  this 
approach  consumes  slightly  more  computer  run  time  than  a  full  step-by-step  calculation 
procedure  it  was  selected  for  two  reasons: 

I  It  allows  easier  alterations  to  each  separate  process  should  changes  in  any 
of  the  local  analyses  be  desired,  e.g.,  a  revised  reaction  kinetics  model.  This 
reduces  the  chances  for  errors  due  to  unforeseen  interactions  between 
sections  of  the  program. 

2.  The  iterative  steps  required  between  the  wake  analysis  and  the  surface 
vaporization  analysis  is  facilitated  if  they  are  isolated  in  the  program. 

The  required  input  consists  of  a  description  of  the  geometry  of  the  streamtube: 

Blockage  ratio 
Flameholder  width 
Flameholder  apex  angle 
Spravbar  to  flameholder  spacing 
Flameholder  to  exhaust  spacing. 
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Those  parameters  which  describe  the  operating  conditions  are: 


Inlet  pressure 
Inlet  velocity 
Inlet  temperature 
Inlet  turbulence  intensity 
Fuel-air  ratio 
Spraybar  fuel  pressure 
Spraybar  fuel  temperature 
Fuel  type 

Vitiation  fuel-air  ratio. 

These  inputs  define  the  required  flamework  for  the  analysis.  The  program  organization  is 
shown  in  Figure  14.  The  various  subroutines  perform  the  analyses  as  follows: 

IN-JKCT  This  performs  the  spray  formation  and  throttling  process  vapor¬ 
ization  analyses. 

ACCEL  This  evaluates  *he  forced  vaporization  and  acceleration  of  the 
droplets  bet  veen  the  spraybar  and  flameholder.  Done  once  per 
droplet  size  'roup. 

COLLECT  This  evaluates  the  collective  rate  of  the  liquid  droplets  onto  the 
surface  of  the  flameholder.  Performed  once  per  droplet  size  group. 

This  evaluates  the  size  of  *he  recirculation  zone  wake  behind  the 
flameholder  and  the  rate  of  entry  of  gaseous  components  into  this 
volume. 

This  evaluates  the  rate  of  vaporization  of  the  liquid  film  from  the 
surface  of  the  flameholder. 

For  the  results  of  RECIRC  and  B3,  this  solves  for  the  wake 
reaction  efficiency  and  temperature. 

This  performs  the  solution  for  the  transverse  penetration  of  the 
turbulent  flame  from  the  wake  shear  layer  into  the  free-stream 
mixture. 

Currently,  the  input  and  output  reflect  the  requirements  to  analyze  a  single  case 
streamline  analysis.  The  output  defines  the  intermediate  results  of  the  various  subroutines  as 
well  as  the  heat  release  profile  between  the  flameholder  and  the  exhaust  nozzle.  The  program 
is  capable  of  multiple  case  execution  limited  only  by  run  time.  This  capability  is  utilized  for 
multiple  streamtube  analysis  and  parametric  studies.  At  present,  there  is  no  graphics  output 
capability,  but  it  could  be  added. 

c.  A.n*lvaet  and  Resulta 

The  combustion  model  as  formulated  for  the  rumble  mechanism  essentially  consists  of 
two  parts.  These  are  a  compositional  analysis  and  a  reaction  analysis.  The  compositional 
analysis  defines  the  manner  in  which  a  combustible  fuel-air  mixture  is  generated  from  the 
injected  liquid  fuel.  The  reaction  analysis  defines  the  manner  in  which  this  mixture  actually 
burns  to  produce  a  final  level  of  combustion  efficiency. 


RECIRC 

m 
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The  analysis  for  the  compositional  portion  is  outlined  below  for  the  basic  case  of  a 
two-dimensional  duct  with  a  bluff  body  stabilizer  and  liquid  fuel  injection  source.  For  ease  of 
development  of  the  solutions,  the  following  assumptions  were  used: 

•  The  injected  liquid  fuel  forms  a  homogeneously  dispersed  spray  which  fills 
the  full  cross  section  of  the  duct. 

•  The  approach  flow  field  is  uniform  in  the  transverse  direction,  i.e.,  trans¬ 
verse  uniformity  of  velocity,  pressure,  and  temperature. 

Under  these  assumptions,  the  form  of  the  compositional  analysis  was  developed.  The 
schematic  of  the  processes  is  shown  in  Figure  15. 


figure  15.  Schematic  of  Wake  Mixture  Formation 


A  certain  portion  of  the  injected  fuel  flowrate  experiences  flash  vaporization  during  the 
injection  process  and  forms  a  uniform  flowrate  of  vapor  phase  fuel.  This  fuel  flowrate  is: 

wtv  =  diTwt.  (14) 

The  remaining  portion  of  fuel  goes  into  a  liquid  droplet  spray  which  undergoes  spray 
vaporization.  The  vapor  flowrate  produced  by  this  process  is: 

w,  =  diK(l~diT)wf.  (15) 

Combining  these  processes  into  one  influence  coefficient  we  may  write: 

wlv  -  d,wf.  (16) 
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We  continue  to  utilize  the  assumption  that  the  vaporized  fuel  is  homogeneously  spread  as  it  is 
produced.  The  next  processes  which  are  evaluated  are  those  concerned  with  the  formation  of 
the  liquid  film  vaporization  rate. 

At  the  flameholder  plane,  the  mixture  consists  of  a  flowrate  of  air  (wa),  vapor  fuel  (wt  ) 


and  liquid  fuel  (w,v)  defined  as: 

wa  =  A  (17) 

WfT  =  0wa  (18) 

Wfv  =  d!WfT  =  0,6  wa  (19) 

Wft  =  w,T  -  Wfv  =  (1-/^1  )(9  wa.  (20) 


This  leaves  (1-di)  liquid  fuel  percentage  available.  Flowing  through  the  flameholder  area, 
then  we  have: 


war  =  Fwa  =  I>  VA 
Wft,.  =  (1-di)  0  waF 
wfvl.  =  0,0  war. 

The  liquid  and  vapor  fuel-air  ratios  at  this  point  are: 


(21) 

(22) 

(23) 


(f/a)v. 


war 

(f/a),  =  (1  -0,)O. 


£i«_war 

w„r 


=  0,8 


(24) 


(25) 


A  portion  of  the  liquid  fuel  is  now  collected  and  vaporized  from  the  surface.  The 
collection  rate  of  liquid  fuel  is: 


w(-  =  02Wfsr  =  (Ml-di)  0  war; 


(26) 


of  this  amount,  /i:l  percentage  vaporizes  from  the  surface  and  recirculates: 

w  =  =  02  0:>(  0i  >  8  waF.  (27) 

The  prevaporized  fuel  recirculates  at  a  rate  of  Kj  percent  which  passes  the  blocked  area: 
w  =  wfvi.  K|  =  di«wa  T  K|.  (28) 

Thus  the  total  recirculated  vapor  fuel  flowrate  is: 

WU  =  <*i«WarK,  +  (l-di)d2d.iflwar  (29) 
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The  air  recirculation  rate  is  K;  of  the  flowrate  through  the  blocked  area: 


w;,w  =  wJ'K,.  (30) 

From  Kquations  16  and  17,  the  wake  fuel-air  ratio  in  vapor  phase  is: 

w, 

</>  -  (31) 

Waw 


,/)  -  di«  +  (1  d.)»  — —  (32) 

Ki 

0  dyds 

—  =  fit  +  (1  -d,)-^—  (33) 

Equation  33  relates  the  wake  vapor  phase  fuel-air  ratio  to  the  overall  fuel-air  ratio  and, 
through  the  influence  coefficients,  to  the  geometric  and  aerothermodynamic  situation  under 
scrutiny.  Each  of  these  may  now  be  evaluated  from  the  available  body  of  combustion  literature 
to  determine  the  overall  form  of  the  compositional  model. 

A  comment  on  the  implication  of  Equation  33  is  in  order  here.  This  equation  determines 
the  mixture  in  the  wake  which  then  undergoes  reaction  to  provide  the  ignition  source  for  the 
mainstream  reaction.  As  with  any  stirred  chemical  reactor,  the  most  efficient  process  for  a 
fixed  volume  occurs  near  the  stoichiometric  fuel-air  ratio.  Since  the  degree  of  this  reaction 
efficiency  strongly  influences  the  overall  augmentor  efficiency,  the  wake  fuel-air  ratio  should 
be  kept  near  stoichiometric. 

If  we  introduce  the  following  values: 

rt,  =  0.20 

li-,  =  0.75 

fa  =  0.50 

K,  =  0.25 

we  find  that  Equation  33  yields: 

</>  (0.75)10.50) 

T  "  °-20  +  <0-8'  025 

140 

0 

Thus,  due  to  the  relative  fuel  concentrating  effect  of  the  flameholder.  the  wake  is  40%  richer 
than  the  overall  system.  If  the  overall  fuel-air  ratio  is  0.050,  the  wake  is  over  stoichiometric  at 
0.070.  Any  further  increase  in  fuel  flowrate  results  in  a  drastic  decrease  in  the  wake  reaction 
efficiency  due  to  the  increase  in  wake  f/a. 

The  second  portion  of  the  analyses  is  the  wake  reaction  efficiency  and  turbulent  flame 
spreading.  These  analyses  will  be  developed  in  the  following  sections. 


(1)  Fuel  Vaporization  Before  Flameholder 


As  developed  in  the  preceding  analysis,  the  percentage  of  the  liquid  hydrocarbon  fuel 
which  vaporizes  prior  to  reaching  the  flameholder  consists  of  two  mechanisms: 

The  formation  of  vapor  due  to  the  throttling  process  of  injection  (di,l- 

The  evaporation  of  fuel  droplets  due  to  forced  convection  (d(  ). 

Thus  we  define  ti,  as: 

d,  d,,  -  (1  d, ,  )^,K.  CM) 

The  process  of  Hash  vaporization  during  the  process  of  injection  is  treated  as  an  adiabatic 
flow  from  a  region  of  high  pressure  and  moderate  temperature  (within  the  sprayring)  to  a 
low-pressure  area  (in  the  augmentor).  'Phis  How  situation  is  analogous  to  the  expansion  valve 
process  and  is  evaluated  from  knowledge  of  the  final  properties  and  operating  conditions. 

The  process  is  evaluated  assuming  adiabatic  expansion  from  known  levels  of  fuel  pressu-e 
and  temperature  within  the  sprayring  to  a  known  pressure  level  in  the  augmentor.  If  the 
sprayring  pressure  is  sufficient  to  maintain  the  liquid  fuel  as  a  saturated  liquid,  the  enthalpy 
level  is  defined  for  the  given  fuel  type  as  a  function  of  the  sprayring  fuel  temperature  only. 


For  adiabatic  injection,  the  final  mixture  enthalpy  equals  the  saturated  liquid  enthalpy  at 
To  The  mixture  enthalpy  and  a  known  static  pressure  will  define  the  quality  of  the  injected 
fuel.  'Phis  process  is  readily  evaluated  from  the  enthalpy  diagram  of  the  particular  iuel  type. 


An  expanded  diagram  of  the  adiabatic  process  for  JIM  is  shown  in  Figure  16  This 
example  is  for  injection  of  250°  F  fuel  into  a  7.5  psia  environment.  The  parametric  results  for 
JF  4  fuel  are  shown  in  Figure  17. 


Fuel  Injection  Temp  -  °F 
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Fif’urr  16.  Constant  Enthalpy  Furl  Injection  of  J P-4  Fuel 


28 


The  remaining  liquid  fuel  is  partitioned  into  the  droplet  size  groups  on  a  5  or  10  group 
equal  flowrate  basis.  For  this,  the  droplet  size  distribution  must  be  known.  Due  to  the  large 
number  of  spray  systems  available  to  the  augmentor  designer,  no  attempt  was  made  to  model 
this  droplet  formation  process.  Rather,  a  distribution  function  is  built  into  the  program  which 
describes  the  size  distribution  as  a  function  of  the  fuel  pressure  drop.  The  curve  used  in  the 
program  represents  the  droplet  distributions  for  variable  area  pintle  spraybars.  Figure  18 
shows  this  function. 

If  5  size  groups  are  to  be  used,  each  one  represents  20%  of  the  liquid  flowrate.  The  sizes 
used  in  the  analysis  would  be  the  10,  30.  50.  70  and  90%  diameters.  Once  these  sizes  are 
known,  the  program  performs  the  solution  to  the  forced  convection  droplet  vaporization 
between  the  spraybar  and  the  flameholder. 

The  basis  for  the  droplet  vaporization  solution  is  the  form  of  the  Nusselt  number  devised 
for  spherical  droplets  by  Ranz  and  Marshall  and  subsequently  improved  by  Preim  and 
Heidmann.  The  assumed  forms  are: 

NuH  2  4  0.6  Re  Pr  (35) 


Nuv. 


2  e  0.6  Re  Sc 


(36) 


tor  heal  and  mass  transfer.  For  the  assumption  of  Lewis  No.  ■=  1,  i.e.,  Pr  -  Sc,  tiiese  become 
identical  The  fuel  droplet  mass  efflux  rate  is  calculated  from: 
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Figure  18.  Typical  Spray  Distribution 


During  the  transient  droplet  heating  period  between  injection  at  Tf  -  T;  and  achieve¬ 
ment  of  the  wet-bulb  temperature,  the  droplet  temperature  is  evaluated  from  the  heat  input 
and  the  diffusion  driven  mass  efflux.  The  difference  between  the  heat  required  to  generate  the 
mass  efflux  and  the  heat  which  actually  reaches  the  surface  of  the  droplet  is  assumed  to  go 
towards  alteration  of  the  droplet  bulk  temperature.  This  assumption  is  essentially  a  statement 
that  the  droplet  internal  circulation  is  sufficiently  rapid  compared  to  the  thermal  input  that 
significant  droplet  radial  thermal  gradients  do  not  exist. 

In  evaluation  of  the  net  heat  flux  into  the  liquid  surface,  the  mass  efflux  blocking  term 
introduced  hv  Preim  is  used.  This  term  evaluates  the  loss  in  net  flux  due  to  the  vapor 
superheating  which  occurs  as  the  evolved  fuel  vapor  achieves  thermal  equilibrium  at  the 
free-strearn  gas  temperature.  The  transient  solution  proceeds  as  follows,  after  evaluation  of  the 
fuel  mass  efflux  rate,  w,  from  Equation  37. 


The  thermal  film  coefficient,  h|,  is: 


b| 


k  Nu(| 
dv 


(39) 


where  Nu(|  is  from  Equation  35. 
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The  net  heat  flux  to  the  droplet  liquid  surface,  q,  is  calculated  with  allowance  for  the 
thermal  blocking  due  to  vaporizing  fuel  heating  (fi): 


q  h,As(Ta  Tt)  ft 


(40) 


where 


d  - 


z 

ez  1 


(41) 


and 

z  ~  Cpv  w/irkdt  Nuh.  (42) 

The  net  amount  of  heat  which  is  available  for  sensible  heating  of  the  droplet  is  the  net 
liquid  surface  flux  minus  the  latent  heat  required  to  generate  the  vapor  mass  flux: 

Aq  ~  q  -  wvA.  (411) 


This  net  flux  is  used  to  raise  the  liquid  bulk  temperature,  T(.  as: 
d'l'v  Aq  Aq 

- - A — -. - -  (44, 

dt  in;  Cpi  4/2  ir  (di/2)  />v  Cpt 

This  change  in  T(  may  be  positive  or  negative  depending  on  the  relative  values  of  w  and  q. 

The  wet  bulb  temperature  is  defined  as  that  liquid  droplet  temperature  where  the 
thermal  net  input  is  just  sufficient  to  generate  the  mass  efflux  at  that  temperature  and  vapor 
pressure.  As  such,  it  is  a  function  of  effective  Nusselt  number  and  local  static  pressure  and 
temperature.  Coincident  with  this  transient  heating  and  vaporization,  the  liquid  droplet  is 
being  accelerated  by  the  faster  free  stream  gas  velocity.  The  acceleration  is  evaluated  from  the 
standard  equation: 


dV( 

HT' 


Pa 

Pi 


<va-v(,)- 


(45) 


Since  the  Reynolds  number  used  in  the  Nusselt  number  formulation  is  defined  as: 


Re 


pA  <Va~Vt) 
H 


(46) 


a  simultaneous  solution  of  Equations  44  to  46  is  required.  In  the  computer  analysis,  this  is 
accomplished  by  a  finite  difference  solution  utilizing  small  time  intervals.  The  ordinary 
differential  equations  are  rewritten  as  delta  terms,  e.g., 

dVc 

dt 


becomes 


At 
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For  each  time  increment,  the  initial  values  are  used  to  calculate  Re,  NuH,  Num,  1\,  etc. 
Equations  37  to  44  are  solved  which  yields  Awv  and  ATp.  Equation  45  is  solved  to  yield  AVt. 
At  the  end  of  this  step,  w,  is  incremented  hy  Awv,  Tt  by  AT(.  and  Vj,  hy  AV(-.  The  average  value 
of  Vv  over  this  At  yields  a  delta  axial  travel  distance. 

The  procedure  is  repeated  in  small  time  steps  until  either  the  total  axial  travel  exceeds 
the  spraybar-to-flameholder  separation  distance  or  the  liquid  is  fully  vaporized.  At  the  start  of 
each  new  time  increment,  the  input  values  reset  to  reflect  the  effect  of  the  previous  step;  e.g., 
/iv  is  evaluated  at  T\  +  A'I\,  dr  is  evaluated  to  reflect  less  liquid  mass,  etc. 

Some  typical  results  for  the  vaporization  portion  of  the  analysis  are  presented  in  Figure 
19.  Figure  19  also  shows  the  transient  hulk  liquid  temperature  for  a  100  micron  droplet  in  a 
400° F,  250  ft/sec  gas  stream.  The  decrease  in  temperature  from  the  initial  80° F  injection 
temperature  for  the  0.5  atmosphere  pressure  reflects  the  fact  that  the  initial  mass  efflux  by 
convection  exceeds  the  available  net  heat  input. 

For  a  situation  described  as  follows: 

V,  =  250  fps 

T„  =  275°F 

T(>  -  80°  F 

P„  =  10  psia 

d|'  -  25  to  125  microns. 

The  amount  of  JP-4  which  will  vaporize  in  a  typical  distance  of  6  inches  is  shown  in 
Figure  20.  These  results  are  typical  for  fan  duct  conditions  at  high  altitude  and  low  flight 
Mach  number. 


Figure  19  Jl’-4  Drnplet  Transient  Temperature 
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Figure  20.  Vaporization  of  JP-i  Fuel  Droplets  vs  Axial  Length 


(2)  Flameholder  Fuel  Collection  Efficiency 

The  output  from  the  ACCEL  subroutine  defines  a  spray  of  droplet  diameters  at  a  plane 
with  the  flameholder  leading  edge.  The  diameters  are  defined  from  the  final  step  in  the  forced 
vaporization  analysis  as  is  the  droplet  axial  velocity.  The  spray  field  is  still  assumed  to  be 
homogeneous  in  the  transverse  direction  in  terms  of  size  and  volume  flowrate  distribution. 


As  the  flow  stream  approaches  the  stabilizer,  the  suspended  droplets  are  unable  to  fully 
follow  the  divergence  of  the  flow  streamlines  and  a  portion  of  them  impinge  on  the  bluff  body. 
The  evaluation  of  this  “capture  rate”  is  done  in  nondimensionalized  terms  as  [}■>: 


W|(.  =  Wfv  1 7L 

(47) 

»  w  Jr 

IF  - 

(48) 

1  wn 

where 

w/(  -  fuel  liquid  capture  rate 

wf(  fuel  liquid  flowrate  at  F/H  plane 
I'  -  F/H  blockage. 

fi>  is  also  the  ratio  of  the  separation  of  the  droplet  capture  limit  streamlines  to  the  F/H 
width: 


wf,  -  in  V,  f>* 

w,i  ;  m  V(  A 

-  M  V,  - 


(49) 
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Substituting  yields: 


h* 


where  6*  ~~  droplet  capture  width.  This  is  shown  schematically  in  Figure  21. 


(50) 


Plane  A 


B2  =  <5*/N  where:  N  =  Flameholder  Width 
2  5*  =  Limit  Streamline 
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Figure  21.  Droplet  Capture  Schematic 


This  analysis  will  be  performed  once  for  each  size  group  exiting  the  vaporization  analysis. 
If  10  groups  are  used  initially,  each  size  represents  10%  of  the  liquid  fuel  at  plane  A,  since  the 
smaller  sizes  will  have  higher  vaporization  rates  between  the  injection  plane  and  plane  A.  The 
overall  collection  rate  is  thus: 
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m, ,  A/i,|. 


(51) 


Evaluation  of  ii>,  requires,  therefore,  evaluation  of  f>*  for  any  value  of  droplet  diameter. 
Additionally,  the  streamline  divergence  rate  will  be  a  function  of  the  flameholder  geometry 
lapex  angle,  blockage).  The  droplet  acceleration  rate  in  the  x  and  y  directions  wid  be  a 
function  of  the  flow  conditions  (pressure,  temperature,  velocity)  and  initial  fuel  droplet 
velocity  (output  from  di,  subroutine). 

An  earlier  analysis  for  collection  into  a  cylindrical  bar  was  performed  by  Langmuir.  The 
collection  efficiency  was  found  to  be  dominated  by  air  velocity  and  droplet  diameter.  The 
influence  of  diameter  is  exceptionally  strong.  Typical  results  are  shown  for  vee-gutter 
flameholders  in  Table  1.  Going  from  25  to  75  micron  droplets  doubles  the  collection  rate  from 
'-:«)%  to  —70%. 


The  effect  reinforces  the  dependence  of  the  low  air  temperature  stabilization  process  on 
the  surface  film  analysis.  As  the  air  temperature  is  reduced,  the  level  of  droplet  vaporization 
is  also  reduced  and  the  droplet  diameter  at  the  flameholder  is  increased.  This  not  only 
increases  the  liquid  flowrate  available  for  collection  but  also  increases  the  percentage  collected. 
This  doubling  effect  produces  a  much  greater  dependence  of  the  wake  fuel-air  ratio  on  the 
surface  film. 


TABLK  I.  (’OI.LKCTK)N  EFFK'IKNCY  vs 
DROPLET  DIAMETER 
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(3)  Fuel  Film  Surface  Vaporization  Evaluation 

The  liquid  droplets  which  contact  the  surface  of  the  flameholder  are  assumed  to  generate 
a  uniform  liquid  film.  This  liquid  film  experiences  heat  transfer  from  the  hot  wake  through  the 
flameholder  as  shown  in  Figure  22. 
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Figure  22.  Flameholder  Surface  Vaporization  Schematic 


The  liquid  film  vaporization  rate  is  used  to  calculate  the  surface  vaporization  efficiency, 
ft  a.  This  efficiency  is  the  percentage  of  the  collected  liquid  fuel  which  is  vaporized  from  the 
flameholder  surface  and  recirculated  into  the  wake  reaction. 


The  controlling  parameter  becomes  a  forced  diffusion  process.  The  controlling  equation 
may  thus  he  written  as: 
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In  these  equations,  the  collected  mass  flowrate,  wt  ,  is  the  result  of  the  ti subroutine.  We 
assume  that  the  liquid  film  total  mass  does  not  change  with  time  unless  the  entry  flowrate 
alters,  i.e.. 


dm 


The  various  f>:m  properties  are  functions  of  the  fuel  type  and  temperature.  Suitable 
curves  are  used  in  th;  program  for  -IFM  and  JP-5  fuel. 

This  situation  is  a  forced  mass  transfer  from  the  surface  represented  by  a  Nusse't 
number  function  driven  by  the  film  vapor  pressure  and  Nusselt  number  form  of  heat  flux  from 
the  wake  into  the  film. 

The  formulation  of  the  rate  of  liquid  fuel  vaporization  from  the  surface  of  the 
flameholder  has  been  programmed  as  a  finite  element  solution  to  the  nonuniform  forced 
diffusion  process.  The  diffusion  process  is  evaluated  from  the  surface  Nusselt  number  for  mass 
transfer: 


wv 


C|  A^psln 


V  IK  -  Pv > 


(58) 


Nu,,,  Dv  MW 
R  Ax  T„ 


(59) 


The  Nusselt  number  is  functionally  identical  for  heat  flux  and  mass  flux  when  the 
Prandtl  number  in  the  thermal  Nusselt  number  is  replaced  with  the  Schmidt  number  for  mass 
transfer. 


The  heat  flux  from  the  recirculation  zone  is  evaluated  directly  from  a  Nusselt  number 
formulation  for  recirculation  wakes  behind  bluff  bodies  as: 


Nuw 


^-=  0.99  Re"5  Pr0:,:'. 


(60) 
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For  a  known  value  of  heat  flux  through  the  flameholder,  the  rate  of  surface  vaporization 
is  evaluated  by  breaking  the  liquid  film  into  10  subgroups.  In  each  subgroup,  the  vaporization 
process  is  assumed  to  be  controlled  by  the  vapor  pressure  corresponding  to  the  mean  liquid 
film  temperature  within  that  subgroup  (see  Figure  23). 

The  heat  flux  into  each  subgroup  is  responsible  for  two  physical  effects.  The  heat  is  used 
to  provide  the  latent  heat  of  vaporization  and  provide  sensible  heat  to  elevate  the  bulk  liquid 
temperature. 

For  the  initial  subgroup,  the  entrance  mass  flowrate  is  the  prorated  portion  of  the  fuel 
collection  rate.  Thereafter,  the  entry  flowrate  consists  of  the  prorated  collection  rate  plus  the 
unvaporized  liquid  fuel  from  the  preceding  subgroup. 

The  entrance  fuel  temperature  for  the  initial  subgroup  is  the  droplet  bulk  temperature  as 
evaluated  in  the  spray  vaporization  subroutine.  For  succeeding  subgroups,  the  initial  fuel 
temperature  is  evaluated  from  the  mixture  of  captured  fuel  at  droplet  temperature  and  liquid 
fuel  film  at  elevated  temperature. 

'[’he  solution  in  each  subgroup  for  the  liquid  temperature  rise  and  vaporization  rate 
requires  an  iteration  process.  The  initial  guess  on  vaporization  is  evaluated  at  the  fuel  entry 
temperature.  The  latent  heat  required  to  accommodate  this  vaporization  is  subtracted  from 
the  heat  flux  and  the  excess  used  to  increase  the  bulk  film  temperature.  A  new  vaporization 
rate  is  calculated  at  a  mean  fuel  temperature  between  entry  and  exit.  The  solutions  for 
vaporization  and  temperature  rise  are  repeated  until  convergence  occurs. 

7'his  process  of  finite  element  solution,  if  the  subdivisions  are  fine  enough,  allows  a 
relatively  simple  ''valuation  of  the  vaporization  rate  for  nonadiabatic  walls  and  a  nonuniform 
liquid  film  temperature.  It  may  be  made  as  exact  as  desired  by  decreasing  the  subgroup  size. 
This  may  be  very  desirable  for  solutions  to  multicomponent  fuels  (such  as  JP-4)  where  the 
latent  heat  of  vaporization  and  surface  vapor  pressure  are  strong  functions  of  the  bulk  liquid 
temperature. 

The  partitioning  of  the  heat  flux,  Equation  56,  and  the  forced  diffusion  process. 
Equations  58  and  59,  may  be  combined  to  yield: 

•  •  _  ,  Num  D,  MW,  ,  p,  , 

q  =  w,;  C,>  (T,  f-Tf  ,)  +  5  (  RAxTa  )  AsPsln  (61) 


where  ('p,  X  and  p^  are  all  functions  of  Tp. 

The  functional  unknown  in  this  group  is  Tp,  the  average  film  temperature.  Once  Tp  is 
known,  the  mass  efflux  is  known  from  Equation  58.  The  form  of  the  efflux  response  to  Tp  is 
such  that  a  critical  film  temperature  exists,  when  pv  =  ps,  where  infinite  flux  exists.  For  JP-4 
and  10  psia  static  pressure  this  is  180°F.  A  very  careful  analysis  is  required  to  iterate  to  a 
successful  answer.  A  typical  result  is  shown  in  Figure  24  for  JP-4  and  flow  conditions  of 
250  fps,  275°F,  10  psia. 
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Figure  24.  Surface  Vaporization  vs  Fuel  Collection  Rate  for  Two  Levels  of 
Wake  Temperature 

(4)  Flameholder  Wake  Recirculation  Ratio 

In  a  manner  analogous  to  the  definition  of  the  collection  efficiency,  we  will  define  a 
recirculation  efficiency,  K(.  This  efficiency  relates  the  mass  flowrate  into  the  wake  of  the 
flameholder  to  the  mass  which  flows  through  the  projected  blocked  area  of  the  flameholder: 

war  =  PaVaNK|.  (62) 

This  recirculation  flowrate  is  the  mass  which  is  transferred  by  turbulent  diffusion  across 
the  free  shear  layers  a  ft  of  the  flameholder.  These  layers  form  the  boundaries  of  the  backflow 
wake  aft  of  the  bluff  body.  If  we  look  at  this  wake  volume  as  a  homogeneous  region  with  mass 
transfer  across  the  boundaries,  the  recirculation  rate  may  also  be  written  as: 


P„V„ 


(63) 


If  the  wake  volume  and  residence  time  may  be  evaluated  as  a  function  of  the  geometric 
and  flow  variables,  then  the  recirculation  rate  may  be  evaluated  from  the  known  variables 
without  resorting  to  the  much  more  difficult  solution  to  the  effective  transport  across  the  shear 
layers. 

The  approach  relates  the  wake  volume  and  residence  time  to  the  controlling  aerodynamic 
and  geometric  variables,  such  as: 

Blockage  ratio 
Vee-gutter  apex  angle 
Flow  Mach  number 
Pressure,  velocity,  temperature. 
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Once  these  are  known,  the  recirculation  rate  calculation  proceeds  as  follows: 


w,r  =  Pa  V,A 

(64) 

V„  =  CV(L/D)(B/D)N-. 

(65) 

If  we  nondimensionalize  the  residence  time  with  respect  to  velocity  and  characteristic 
dimension,  flameholder  width,  we  have: 


(66) 
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Now,  from  Equations  65  and  67,  we  have: 

,,aVflCv(L/D)(B/D)N 
...  —  ~ 


(68) 


From  this  and  Equation  62,  we  may  write: 

K,  -  Cv(L/I))(B/D)(t  )  (69) 

The  value  of  the  shape  coefficient,  Cv,  was  determined  by  assuming  that  the  wake  was 
similar  to  a  two-dimensional  ellipse  which  is  truncated  by  the  flameholder  at  its  forward  edge. 
The  value  used  is  0.80. 

The  data  presented  were  reduced  to  a  series  of  curves  which  describe  the  effect  of  various 
parameters  on  recirculation  rate.  The  functional  groupings  and  the  corresponding  figures  are: 


L/D 

vs  Blockage 

(Figure  25) 

L/D 

vs  Apex  Angle 

(Figure  26) 

B/D 

vs  Blockage 

(Figure  27) 

B/D 

vs  Apex  Angle 

(Figure  28) 

t' 

vs  Blockage 

(Figure  29) 
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Figure  25.  Wake  L/D  vs  Hlockagt •  Ratio 


r  =  Blockage  Ratio  -  D'Less 
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Figure  29.  Wilke  Residence  Time  ixs  Blockade  Ratio 


(5)  Recirculation  Zone  Reaction  Kinetics 

The  recirculation  zone  wake  of  the  bluff  body  stabilizer  is  analyzed  by  assuming  that  it 
behaves  similarly  to  a  well-stirred  reactor  with  the  volume  and  mass  entry  rate  known  from 
the  results  of  the  REC1RC  analyses. 


Analyses  of  well-stirred  reactors  have  been  presented  by  numerous  investigators  and  the 
analogy  to  actual  combustors  and  bluff  body  wakes  studied.  The  basic  thesis  of  these  studies 
is  that  the  performance  of  the  reactor  may  be  evaluated  from  a  balance  between  the  mass 
entry  rate  and  the  kinetic  conversion  rate.  For  the  purpose  of  these  studies,  we  assume  that 
the  entire  wake  volume  is  available  for  reaction  and  that  the  mixing  is  very  rapid.  These 
assumptions  result  in  an  optimistic  evaluation  of  the  kinetic  limits  which  we  will  correct  later. 

The  following  analysis  follows  the  development  of  Kretschmer  and  Odgers  for  lean  wake 
fuel-air  ratios.  The  general  form  of  the  reaction  of  hydrocarbon  fuels  proceed1'  in  essentially 
two  major  steps.  The  first  is  the  pyrolysis  and  partial  oxidation  of  the  virgin  fuel  to  form 
short  lived  intermediate  species.  At  the  end  of  this  stage,  the  maximum  concentration  of 
carbon  monoxide  is  present.  The  second  stage  is  the  oxidation  of  carbon  monoxide  to  form 
carbon  dioxide.  This  latter  step  is  much  slower  and  serves  as  the  rate-controlling  process.  The 
slower  CO  oxidation  results  in  the  very  close  similarity  of  reaction  rates  for  a  wide  range  of 
saturated  hydrocarbons. 


The  general  form  of  the  conversion  equation  max  be  written  as: 


dm,, 

dt 
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for  gaseous  flow  into  a  well-stirred  reactor,  this  becomes: 
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(71) 


For  the  assumed  single-step  reaction  process  postulated  here,  the  reaction  mass  balance  is  (for 
propane  fuel): 


(pt'iiHfi  +  50-j  l  5mN-_i  — .  3y<  CO-,  t  4v<  M.O  t  (72) 

(</>— y< )  +  5(1  y<)  02  +  OmN,. 

Also,  a  linear  efficiency  vs  temperature  function  is  assumed: 

T  =  Tj  +  <AT.  t73) 

From  these  equations,  the  stirred  reactor  loading  capability  may  be  written  as: 

A  k(m  +  l)  1 5(1 ->■»)“  (./.  -y< )"  a  e  ('^1i^<AT) 

V„p"  =  R,ny<  |5(m+ 1 )  +  ./.  +  y<|"  |T,+:ATr  "  '  ('4) 

Based  on  comparison  of  predicted  results  with  available  stirred  reactor  data,  Reference  21 
recommends  the  following  values  for  this  reaction: 

n:  for  <l>  <  1 .  n  -  2</> 

for  (/>>l,n=2A/> 
a:  a--n/2. 

This  yields: 

A  1.29  x  101"  tm-tl)  [5(  1  -y<|';’  [fi-yi)*  e  (VlTi+tAT) 
pp2'*  ~  W.08206)-:’  y<l5(m+l )  +  </>  +  yep  |T,  +.ATJ 

lor  lean  mixtures. 


While  the  results  presented  above  are  from  one  reference,  similar  results  have  been 
obtained  by  others  for  the  same  problem. 

Longwell  and  Weiss  (Reference  4)  present  the  following  results  for  lean  and  rich 
operation. 


Lean: 


V„p,H 

Rich: 

A 
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(77) 


where 


k,  =  1.67  X  10"’  litre0H/°Ku 5  gm-mole0H  sec 
k,  =1.11  X  10"  litre0 l7°K"r’  gm-mole0K  sec 
E  =  42,000  gm-cal/gm-mole. 

The  solution  procedure  utilized  in  the  WAKE  subroutine  utilizes  the  results  of  RECIRC 
as  a  definition  of  the  entry  rate,  A,  and  the  zone  Volume.  V„.  The  equivalence  ratio  is  from 
Equation  33  and  the  results  of  d  b  d  2,  d  3  and  RECIRC.  The  ideal  temperature  rise,  AT,  is 
read  from  curves  of  AT  versus  wake  fuel -air  ratio,  inlet  temperature,  and  static  pressure  for 
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this  specified  fuel  type.  An  iterative  solution  procedure  solves  for  the  efficiency,  <,  where  the 
known  value  of  A/V„p‘  agrees  with  the  predicted  value  as  a  function  of  *.  The  complexity  of 
Equation  76  or  77  requires  this  sort  of  reverse  solution. 

The  behavior  of  the  solution  is  shown  in  Figure  30  for  inlet  temperatures  of  300  and 
400°K  and  equivalence  radios  of  0.8  and  1.0.  The  wake  efficiency  decreases  at  an  increasing 
rate  until  the  detuv  slope  becomes  infinite.  At  this  point,  the  wake  reaction  process  is  said  to 
have  blownout.  If  we  plot  the  locus  of  the  blowout  points  as  A/V,,pJ  versus  equivalence  ratio, 
we  obtain  a  classical  blowout  curve.  This  curve,  except  for  magnitude,  represents  the  classical 
results  for  gaseous  fuel  data  from  wake-stabilized  flames.  This  is  shown  in  Figure  31. 
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Figure  3(1.  Stirred  Reactor  Kinetic  Efficiency  vs  Loading  Rate 

This  analysis  tends  to  overestimate  the  limits  of  blowout  velocity  when  compared  to 
available  data.  To  reconcile  this,  comparisons  were  made  between  predicted  limits  and  actual 
limits  for  available  data.  The  correlating  ratio  between  predicted  and  actual  was  3.55,  i.e.: 


3.55. 


(78) 


In  the  analysis,  the  calculated  recirculation  rate  per  unit  volume  is  multiplied  by  this 
factor  before  solving  for  t.  In  this  manner,  the  predicted  limits  and  available  data  are 
numerically  equal,  and  the  efficiency  response  towards  blowout  follows  the  theoretical  curve. 

The  reason  for  this  error  is  most  likely  either  imperfect  mixing  or  lack  of  full  utilization 
ot  the  wake  volume  for  reaction.  The  approach  used  above  will  account  for  either  one  of  these. 
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(6)  Turbulent  Flame  Spreading  Rate 

The  turbulent  flame  propagation  into  the  unreacted  free-stream  is  initiated  in  the  shear 
layers  of  the  wake.  The  model  used  (subroutine  ST)  relates  the  local  turbulent  flame  speed  to 
the  local  aerothermodynamic  conditions  and  performs  a  finite  difference  integration  of  the 
(lame  front  penetration  starting  in  the  wake  and  proceeding  to  the  exhaust  nozzle. 

For  the  purposes  of  current  analysis,  the  following  assumptions  were  made: 

1.  Uniform  air  How  profiles 

2.  Uniform  fuel-air  ratio 

,1.  Incompressible  acceleration  of  free  air  velocity  by  the  flameholder 

blockage  with  no  induced  profile 

4.  Known  wake  size  and  reaction  efficiency 

5.  Two-dimensional  ducted  flame. 

The  schematic  of  the  situation  which  is  analyzed  is  shown  in  Figure  112. 

The  approach  flow,  at  known  levels  of  pressure,  temperature,  velocity  and  fuel-air  ratio, 
is  accelerated  by  the  blockage  of  the  flameholder  to  velocity  U,  where: 


where: 

U  —  Velocity  at  flameholder  tip 

Va  ~  Approach  velocity 

I'  ~  Blockage  ratio. 

At  this  point,  Station  1  of  Figure  31,  an  induced  turbulence  level  is  calculated  from: 

•  [  •  •  (tt)  i]M 

This  equation  ( Reference  12)  relates  the  turbulence  intensity,  to  the  blockage  ratio  and  the 
flameholder  zero  blockage  drag  coefficient,  Cj. 

At  this  location,  the  turbulent  flame  velocity  calculations  are  initiated.  The  equation  used 
for  the  local  flame  speed  is  the  Karlovitz  equation: 

St  =  Sg  f  (2u'SiTJ  (81) 

where: 

St  ~  Turbulent  flame  speed,  ft/sec 

S4  ~  Laminar  flame  speed,  ft/sec 

ij’  ~  RMS  turbulence  velocity,  ft/sec. 

'Che  value  of  u'  is: 


ii  -  U. 
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Air  Velocity  Accelerated 
By  Blockage  -  Uniform 
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Figure  32.  Schematic  of  Flame  Spreading  Anahsi: 


Aii  additional  term  is  required  to  relate  the  resultant  Maine  speed  to  the  level  of 
reel reuiat ion  zone  reaction  efficiency.  As  the  wake  efficiency  and  temperature  decrease,  local 
areas  in  the  flame  sheet  appear  where  ignition  does  not  occur.  Since  the  model  assumes  a 
continuous  Marne  sheet,  we  model  these  local  areas  of  no  ignition  by  reducing  the  overall  Marne 
speed.  A  reduced  Marne  speed  results  in  lower  overall  combustion  efficiency  which  is  similar  to 
the  effect  of  locally  zero  ignition.  The  model  uses  the  following  equation: 

St'  -  St  •  rAv.  (83) 

The  initial  value  for  the  augmentor  efficiency  is  the  wake  reaction  level  on  a  mass 
weighted  basis.  Expressed  as  an  equation  this  is: 

mr 

>7,  =  Vk  '  ~  (84) 

m., 

where: 

>/,.  ~  Initial  efficiency 

?7W  ~  Wake  efficiency 

mr  ~  Wake  mass  Mowrate 

ma  ~~  Total  duct  Mowrate. 

The  type  of  Marne  utilized  in  this  model  is  a  zero  thickness  Marne  which  separates  a  region 
of  unreacted  propellants  from  a  region  of  completely  reacted  products.  From  this  setup  the 
average  local  augmentor  efficiency  is  simply  the  ratio  of  the  transverse  Mame  penetration,  Ay, 
to  the  duct  width,  W  (see  Figure  32). 

To  be  consistent,  the  transverse  location  of  the  Marne  front  at  the  initial  calculation 
station  is  taken  to  be: 

Ay„  =  i/(.  •  W.  (85) 

This  value  is  assigned  to  the  first  axial  station,  which  is  assumed  to  occur  halfway  down 
the  length  of  the  recirculation  zone.  From  visual  observations  of  wake-stabilized  Mames,  this  is 
the  approximate  location  of  transverse  Marne  initiation. 

From  this  location  downstream  to  the  exhaust  nozzle,  the  Mame  front  transverse  location 
is  calculated  by  a  finite  difference  integration  of  the  local  Mame  speed.  Several  axial  profiles 
are  introduced  as  the  integration  proceeds.  These  are: 

•  The  turbulence  intensity  is  decayed  from  the  value  generated  at  the  aft 
Mameholder  lip  (Equation  80)  at  a  rate  inversely  proportional  to  the  square 
root  of  axial  distance  over  an  effective  jet  length.  The  final  value  is  set  at 
the  initial  turbulence  level.  The  effective  jet  length  is  set  at  10  L/D  where 
the  I)  is  the  open  area  distance  between  adjacent  Mameholders. 

•  The  velocity  of  the  unreacted  fuel-air  mixture  is  retained  at  the  level 
generated  at  the  Mameholder  lip.  Measured  profiles  from  several  ducted 
Mame  test  rigs  support  this  assumption. 

•  A  term  is  introduced  which  relates  the  local  Mame  speed  to  the  local 
average  duct  combustion  efficiency,  peaking  at  50%.  This  treats  the 
counteracting  inMuences  of  reduced  heat  loss  as  efficiency  increases  and 
the  reduced  free  oxygen  concentration.  Local  rates  which  follow  roughly  a 
sine  wave  function  have  been  reported  from  duct  data. 
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Analysis  of  the  terms  utilized  for  evaluation  of  the  laminar  flame  speed  term,  ST  has 
resulted  in  the  following: 


Sv 


(</>  ) 


(86) 


where: 

Si  =  laminar  flame  speed  at  1  atm  and  540  deg 
</>  equivalence  ratio 
T„  -  air  temperature  °F 
X„.,  ~  oxygen  mole  fraction. 

'l'he  influence  of  pressure  is  indeterminate  at  this  time  and  has  been  incorporated  as  y  I' 
for  subatmospheric  data  and  no  influence  for  pressures  above  1  atmosphere. 

d.  Phase  I  Results 


The  first  efforts  to  improve  the  combustion  model  were  directed  at  the  solution 
methodology  for  the  fan  stream  flamehoider  wakes.  The  solution  technique  for  a  given  fan 
streamtube  utilizes  two  sets  of  equations  with  the  flamehoider  wake  temperature  and 
flamehoider  wake  fuel-air  ratio  as  independent  variables. 


The  first  equation  set  is  called  the  composition  solution.  For  a  given  geometry  and 
operating  point,  the  wake  fuel-air  ratio  (</>)  is  related  to  the  streamtube  fuel-air  ratio  (H)  by  the 
equation: 


£ 

a 


di  +  <l-di> 


Ihh 

K; 


In  this  equation,  all  terms  are  independent  of  the  wake  temperature  except  d,.  Since  d,;  is 
controlled  by  the  heat  flux  through  the  flamehoider  surface,  the  wake  fuel-air  ratio  will 
increase  as  the  wake  temperature  increases. 

l'he  solution  technique  arbitrarily  increases  the  wake  temperature  through  a  range  from 
10t)t)°F  to  5000°F  and  calculates  the  resultant  value  of  wake  fuel-air  ratio  as  d:;  increases.  Such 
solutions  are  shown  in  Figure  88  for  one  geometry  and  operating  point  for  three  levels  of 
streamtube  fuel-air  ratio. 

The  vertical  portion  of  the  composition  solution  curve  for  0.020  streamtube  fuel-air  ratio 
occurs  because  /|,  reac  hed  a  value  of  1.0  at  about  2600°F.  Since  all  the  surface  fuel  film  was 
vaporized  here,  further  increases  in  wake  temperature  resulted  in  no  further  increase  in  wake 
fuel  entry  or  fuel-air  ratio. 

The  second  equation  set  reverses  the  order  of  fuel-air  ratio  and  temperature  to  produce 
a  curve  of  wake  temperature  versus  fuel-air  ratio.  In  this  case,  the  wake  fuel-air  ratio  is 
arbitrarily  increased  from  a  very  low  value.  The  wake  reaction  is  evaluated  at  this  fuel-air  ratio 
and  the  given  set  of  geometry  and  inlet  conditions.  For  the  initial  low  values  of  fuel-air  ratio, 
no  wake  kinetic  solution  will  exist.  The  first  fuel-air  ratio  value  which  produces  a  wake 
reaction  solution  establishes  the  lean  limit  for  the  given  conditions.  The  solution  technique 
continues  to  increase  the  wake  fuel-air  ratio,  establishing  a  wake  temperature  for  each  fuel-air 
ratio,  until  a  value  of  wake  fuel-air  ratio  is  tried  which  does  not  yield  a  solution.  This 
establishes  the  wake  rich  limit. 
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The  curve  which  results  from  this  approach  is  shown  in  Figure  :{;i  labeled  kinetics 
solution. 


Wake  Fuel-Air  Ratio 


Figure  :i:i.  Duct  Stream  Flameholder  Wake  Solution 

The  final  step  in  the  solution  for  the  fan  duct  streamtube  is  to  locate  the  simultaneous 
solution  of  these  two  procedures.  The  intersection  defines  the  possible  solution  to  the  fan  duct 
streamtube  wake  fuel-air  rat i<  and  wake  temperature.  It  is  also  possible,  however,  that  there  is 
no  intersection  of  these  two  curves.  The  composition  curve  may  miss  the  kinetics  curve  on 
either  the  lean  or  riel,  side.  These  situations  reflect  the  use  of  input  variables  for  which  no 
stable  llame  is  possible.  Since  it  is  not  possible  to  a  priori  determine  the  lean  and  rich  limits 
for  all  configurations,  the  intersection  solution  must  recognize  that  failure  to  intersect  is  not  an 
error  but  a  valid  blowout. 

The  solution  vhich  is  used  in  the  computer  program  was  found  to  have  several  serious 
flaws  which  were  directl\  responsible  for  inconsistent  answers  from  the  program.  The  program 
required  that  variations  in  overall  efficiency  with  changes  in  inlet  pressure,  temperature, 
velocity  and  fuel-air  ratio  should  be  smooth  curves.  The  solution  technique  originally  em¬ 
ployed  resulted  in  a  variation  as  shown  in  Figure  34. 

The  problems  were  traced  to  several  defects.  First,  the  two  sets  of  curves  were  stored  in 
arrays  which  were  not  dimensioned  large  enough  to  store  the  full  range  of  data  at  some  points. 
This  resulted  in  a  true  solution  which  existed  beyond  the  range  of  the  array,  but  was  solved 
using  the  last  array  value.  As  a  result,  the  “solution”  would  vary  depending  on  exactly  where 
the  array  started. 

The  second  problem  was  related  to  the  solution  for  wake  reaction  efficiency  as  a  function 
of  the  wake  aerodynamic  loading  parameter  (V/PsT1  ').  This  function,  for  a  given  loading  and 
fuel-air  ratio,  has  two  possible  solutions,  one  of  which  is  invalid.  Therefore,  provisions  were 
made  to  assure  selection  of  the  proper  solution.  Unfortunately,  the  provisions  covered  all 
possibilities  but  one.  That  possibility  was  when  an  initial  guess  on  the  solution  was  greater 
than  a  failure  solution  but  within  the  convergence  tolerance.  This  problem  was  eliminated  by- 
using  the  slope  of  the  function  at  convergence  to  determine  if  the  solution  is  valid. 


A  reaction  kinetics  solution  curve  was  generated  by  incrementing  fuel-air  ratio  and 
solving  for  wake  temperature  until  a  point  failed.  From  that  point,  the  f/a  was  reduced  in  very 
small  increments  until  a  point  solved.  'I'his  was  defined  as  the  rich  limit.  The  problem  occurred 
when  the  rich  limit  point  was  actually  the  last  solved  point  (i.e.  f/a  reduced  a  full  increment). 

In  this  case,  a  round-off  error  occurred  in  the  seventh  decimal  place  in  f/a  causing  the 
curve  interpolating  routine  to  return  erroneous  temperature  information  when  evaluating  near 
the  rich  limit.  This  problem  was  resolved  by  eliminating  this  extra  solution  information 
whenever  the  error  occurred. 

Also,  the  slight  waviness  in  the  results  was  eliminated  bv  tightening  the  tolerance  on  all 
iterative  solutions  from  0.005  to  0.0001.  Graphic  examples  of  the  results  of  these  efforts  are 
shown  in  Figures  54  through  57.  The  data  points  used  are  from  the  engine  data  taken  under 
contract  F55615-76-G-2024  and  presented  in  AFAPL-TR-78-82.  Variations  are  shown  plotted 
versus  static  pressure  (PJ  and  fuel-air  ratio  (f/a). 


Figure  35.  Revised  Combustion  Model  Results,  Fan  Stream  Case 
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Figure  37  Sample  Case:  Data  Pt  195  B/M 
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The  revised  solution  technique  also  eliminated  the  problem  of  a  blowout  prediction  (/; 

0)  in  between  two  valid  solutions.  Now,  once  the  flame  blows  out  it  stays  out. 

The  next  area  of  investigation  was  to  determine  why  the  surface  vaporization  model 
would  not  execute  for  approach  air  temperatures  in  excess  of  approximately  .‘fr>0°F.  The 
original  surface  vaporization  model  was  formulated  anti  calibrated  for  approach  air 
temperatures  in  the  range  of  200  to  250°F.  One  of  the  heat  transfer  modes  which  was  included 
in  the  model  was  heat  transferred  from  the  airstream  to  the  surface  fuel  film  on  the 
flameholder. 

As  test  cases  w'ere  run  to  evaluate  the  range  of  applicability  of  the  model,  it  was  found 
that  premature  rich  wake  blowouts  were  predicted  above  I150F  and  by  500F  virtually  no  stable 
flame  was  possible.  This  result  is,  of  course,  in  opposition  to  observed  (lame  stability  behavior. 
The  investigation  into  the  reson  for  this  failure  focused  on  the  formulation  for  the  vaporization 
rate  of  a  liquid  exposed  to  high  velocity  airflow  with  subsurface  heating. 

'The  results  of  this  study  showed  that  the  Nusselt  number  formulation  for  the  heat 
transfer  from  the  air  flow  to  the  liquid  fib-’  was  in  error  on  the  high  side.  The  result  of  this 
error  was  that  the  majority  of  the  heat  wi  ,h  was  responsible  lor  the  liquid  vaporization  was 
due  to  the  air  flow  and  not  due  to  the  heat  from  the  flameholder  wake.  As  the  air  temperature 
increased,  the  degree  of  surface  vaporization  increased  (larger  d, I  and  the  wake  fuel-air  ratio 
exceeded  the  rich  limit. 

The  equations  were  revised  to  reduce  the  dependence  on  the  air  heat  transfer  mode.  This 
reduced  the  amount  of  surface  vaporization  considerably  and  resulted  in  converged  solutions 
to  well  beyond  8(M)°F.  However,  another  problem  arose.  File  converged  solutions  were  now 
quite  a  bit  leaner  and  very  large  values  of  approach  fuel-air  ratio  were  required  for  blowout, 
due  to  the  low  percentage  of  heat  flux  through  the  wake 

Additionally,  when  a  fuel  other  than  .JP4  was  used,  specifically  -JF5,  the  model  refused  to 
converge.  The  lower  volatility  of  JP5.  coupled  with  the  lower  heat  flux  into  the  liquid,  did  not 
yield  enough  tporized  fuel  to  achieve  wake  fuel-air  ratios  above  the  lean  limit.  The  situation 
was  thus  that  there  was  apparently  another  mode  of  transfer  of  liquid  fuel  into  the 
flameholder  wake  which  was  not  accounted  for  originally  and  the  excessive  air  heat  flux  gave 
about  the  correct  level  of  extra  fuel  at  about  250° F. 

This  missing  mechanism  was  discovered  in  a  paper  by  (later  and  I/Ecuyer  (Reference  1) 
on  liquid-film  cooling  mass  transfer.  Two  modes  of  liquid  mass  transfer  to  the  local  shear  layer 
are  postulated.  The  first  is  vaporization  as  was  used  in  the  model.  The  second  employed  the 
direct  transfer  of  liquid  from  the  film  to  the  shear  layer  due  to  roughness  of  the  liquid  surface. 

This  mechanism  was  added  to  the  calculation  routine  for  in  the  model.  Since  there  is 
little  dependence  on  the  air  temperature,  excessive  wake  enrichment  is  avoided.  Since  the 
direct  transfer  is  also  not  volatility  dependent,  the  extremely  lean  results  previously  obtained 
with  Jl’ft  fuel  shifted  to  reasonable  wake  fuel-air  ratio  values. 

A  significant  problem  was  discovered  in  the  combustion  model  during  the  course  of 
evaluation  of  multiple  streamtube  configurations.  A  streamtube  was  used  in  the  core  stream  to 
represent  a  location  in  the  flameholder  where  no  reaction  was  present.  To  model  this,  the  user 
input  an  extremely  narrow  flameholder,  which  should  have  been  past  the  blowout  limit.  The 
model  predictions,  however,  resulted  in  stable  flames  to  well  below  experience  limits.  The  wake 
stability  model  was  severely  overpredicting  the  blowout  limits  at  the  higher  levels  of  approach 
gas  temprature  which  exists  in  the  turbine  discharge  air  flow  or  core  region.  A  reexamination 
of  the  formulation  for  the  blowout  limits  was  started. 


The  formulation  for  the  core  region  flameholder  wake  consisted  of  a  recirculation 
calculation  and  the  wake  kinetics  calculation.  The  basic  thesis  is  that  the  wake  region  could  be 
treated  behaviorally  as  a  stirred  reactor  whose  volume  was  equal  to  the  recirculation  zone  size 
and  whose  mass  loading  rate  was  equal  to  the  recirculation  rate.  This  reactor  was  assumed  to 
have  a  fuel-air  ratio  equal  to  the  injected  fuel  flowrate  divided  by  the  approach  air  flowrate. 

The  problems  which  were  discovered  included  ar.  inappropriate  temperature  exponent  on 
the  recirculation  rate  which  reduced  the  effective  reactor  loading  at  elevated  temperatures. 
This  reduced  loading  calculation  permitted  higher  levels  of  wake  reaction  efficiency  and  thus 
improved  blowout  limits.  Additionally,  the  analysis  originally  did  not  account  for  the  effect  of 
the  vitiation  of  the  core  stream  approach  tlow  upon  the  reaction  zone  equivalence  ratio  or 
adiabatic  flame  temperature. 

The  recirculation  rate  and  the  reactor  loading  rate  were  originally  related  as  follows: 

•  The  reactor  loading  parameter,  after  Longwell  (Reference  4)  is: 


A 
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where: 

mass  entry  rate 
reactor  volume 
static  pressure 
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•  The  recirculation  rate  is: 


A  =  K,  ,.VrN  =  K,  -j—-  VK  N 


•  The  reactor  volume  is: 

V,,  =  Constant  ■  (L/d)(B/d)NJ 


Combining  yields: 


Constant  *  K,  VgN  ■  (L/d)(B/d)NJ  ' 
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•  An  air  viscosity  term  was  added  to  account  for  the  assumed  dependence 
of  recirculation  rate  in  approach  viscosity  through  temperature,  yielding 
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This  relationship  has  the  form  of  the  classical  blowout  correlations  and  was  thought  to  be 
correct.  What  was  overlooked  was  the  allowable  limit  on  reactor  loading  would  increase  as  the 
approach  temperature  increased.  Since  blowout  occurs  when  the  calculated  recirculation  or 
loading  exceeds  the  limit  value,  the  result  of  this  oversight  was  to  raise  the  blowout  limits 
quite  rapidly  at  elevated  temperatures.  The  model  would  predict  a  stability  limit  which  would 
correlate  with  Tg  as  opposed  to  the  experimental  data  range  of  TK  12  to  TR  1  ’. 

The  corrections  were  twofold.  First,  the  assumed  viscosity  term  was  deleted.  Secondly,  an 
effect  of  approach  temperature  on  the  wake  volume  was  noted  and  included.  The  addition  of 
these  two  corrections  resulted  in  a  loading  parameter  as: 
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when  coupled  with  the  T1  *  increase  in  the  limit  value  the  stability  correlation  term  became 
V 

pNT  ' 

This  was  better,  but  still  overpredicted  the  stability  limit  increase  due  to  Tg  increase  when 
compared  with  available  data. 


The  final  correction  was  made  when  it  was  noticed  that  the  flameholder  wake  kinetics 
analvsis  was  being  performed  using  a  fuel-air  ratio  which  was  based  on  the  injected  fuel 
flowrate  only  and  ignored  the  ef'ect  of  the  vitiation  which  occurs  in  the  main  combustor.  This 
effect  was  modeled  in  the  wakes  as  an  increase  in  the  diluent  fraction  of  the  recirculated  air 
flow  with  te  reaction  fuel-air  ratio  based  on  the  available  oxygen  mole  fraction. 

'Phis  revised  analysis  had  the  effect  of  offsetting  the  increase  in  stability  limit  as 
temperature  increased  since  the  available  oxygen  decreased.  This  effect  is  noted  only  for  the 
situation  where  the  approach  gas  temperature  is  raised  by  direct  vitiation. 

To  evaluate  me  success  of  this  revised  model,  a  series  of  stability  limit  calculations  were 
performed.  The  peak  of  the  blowout  curve  was  plotted  as  a  function  of  vitated  approach  gas 
temperature  The  resultant  curve  was  then  compared  to  several  gaseous  stability  limits  from 
open  literature  and  internal  data.  It  was  expected  that  the  resultant  curve  would  have  a 
temperature  exponent  near  the  value  of  1.7. 

The  result  is  shown  in  Figure  88.  The  stability  limit  does  not  show  the  expected  trend 
but  neither  did  the  data.  The  model  does  show  the  same  trend  as  the  data  indicated.  The 
conclusion  is  that  the  model  represents  a  realistic  evaluation  of  the  fundamental  physics  of 
bluff  body  flame  stabilization.  Also,  at  least  for  vitiated  air  flow,  there  is  no  single  value  of  the 
exponent  on  temperature. 

The  final  modification  which  was  made  to  the  combustion  model  was  in  the  area  of  the 
initialization  of  the  turbulent  flame  speed  from  the  shear  layers  of  the  recirculation  zone.  This 
interaction  is  necessary  to  couple  the  effects  of  inlet  dynamics  on  wake  conditions  to  the 
ultimate  level  of  overall  efficiency.  The  working  thesis  is  that  as  the  inlet  conditi  n;  are 
altered,  e.g.,  pressure  decrease,  velocity  increase,  the  wake  reaction  efficiency  also  charges. 
This  change  in  inlet  conditions  affects  the  turbulent  flame  speed  directly  through  the  flame 
propagation  model  and  indirectly  through  the  process  of  ignition  from  the  wake. 
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If,  for  example,  the  inlet  velocity  were  increased,  the  wake  efficiency  and  thus  shear  layer 
temperature  would  decrease.  This  would  result  in  a  lower  energy  level  for  ignition  of  the  bulk 
gas  How  and  a  greater  potential  for  ignition  failure.  The  combustion  analysis  models  this  effect 
through  an  equation  which  relates  the  initial  turbulent  flame  speed  to  the  wake  efficiency  and 
shear  layer  temperature. 

The  original  analysis  used  a  term  which  decreased  the  turbulent  flame  speed  by  the 
linear  factor  of  wake  efficiency  and  surface  vaporization  percentage.  The  vaporization  was  used 
to  account  for  the  quenching  effect  of  liquid  fuel  present  in  the  shear  layer.  This  simple  model 
did  not  result  in  the  dramatic  decrease  in  overall  efficiency,  as  the  blowout  was  approached, 
which  was  required  to  drive  the  rumble  instability. 

A  revised  coupling  model  was  used  in  the  final  version  of  the  program.  This  model  uses 
a  strong  nonlinear  functional  relationship  between  shear  layer  temperature  and  initial  flame 
speed.  This  is  felt  to  be  a  more  accurate  description  of  the  initial  steps  in  the  ignition  process. 
The  results  of  this  revised  model  are  shown  in  Figure  39.  The  increased  dependence  of  the 
overall  efficiency  on  the  wake  conditions  is  seen  as  the  stronger  decrease  rate  near  blowout. 
The  ultimate  stability  limit  is  not  altered  by  this  revision,  since  the  wake  blowout  limit  is 
unchanged.  The  rumble  onset  point  is  moved  away  from  the  blowout  point  sufficiently  to  allow 
a  less  critical  evaluation  by  the  user. 
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Figure  119  Sample  Case:  Data  Pt  504  Uhl 
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A  further  modification  in  the  model  was  made  to  calculate  the  wake  recirculation  time 
constant  which  is  used  in  the  stability  analysis.  This  constant  represents  the  time  lag  between 
which  a  dynamic  flow  perturbation  is  seen  at  the  flameholder  trailing  edge  and  the  effect  is 
transferred  to  the  turbulent  flame  propagation  rate. 

The  original  model  assumed  that  the  effects  of  flow  perturbations  were  transferred  across 
this  interface  instantaneously.  The  revised  model  accounts  for  the  time  which  is  required  for 
the  wake  region  to  respond  to  a  flow  shift.  This  time  constant  is  already  calculated  as  the  wake 
residence  time  in  the  combustion  recirculation  rate  solution. 

The  revision  consisted  of  evaluating  the  average  value  of  this  constant  over  the 
streamlines  used  and  placing  the  value  in  storage  for  use  by  the  stability  analysis. 


e.  Computer  Program  Description 

The  computer  code  is  currently  set  to  analyze  a  conventional  turbofan  augmentor  with 
vee-gutter  flameholders.  The  liquid  fuel  is  injected  through  a  spraybar  located  a  specified 
distance  upstream  of  the  flameholder.  The  combustion  process  is  stabilized  by  the  bluff  body 
recirculation  zone,  and  a  turbulent  flame  sheet  propagates  into  the  approach  fuel-air  mixture. 
The  geometry  is  two-dimensional  and  would  represent  one  streamtube  of  a  multiple  stream 
augmentor  system. 

The  various  analyses  which  were  developed  previously  and  referenced  in  the  following 
paragraphs  are  written  as  subroutines  in  the  computer  program.  This  results  in  a  modular 
program  with  a  supervisory  MAIN  executing  subroutines  as  required.  This  approach  allows 
easy  modification  of  the  various  analyses  without  disturbing  the  overall  program  operational 
logic. 


The  program  first  reads  the  input  in  NAMELIST  format.  From  the  input  the 
two-dimensional  model  format  is  set  up.  Total  air  and  fuel  flowrates  are  calculated  and  the 
core  stream  or  duct  stream  option  is  exercised.  The  following  description  details  the  duct 
stream  analysis  procedure. 

From  the  fuel  conditions  in  the  spraybar  and  the  flow  field  conditions  in  the  duct,  the 
degree  of  flash  vaporization  which  occurs  during  injection  is  calculated.  This  percentage  of  fuel 
is  allocated  to  the  initial  vapor  phase. 

The  fuel  which  remains  liquid  is  assigned  to  five  equal  mass  flowrate  groups.  These 
groups  are  each  assigned  a  mean  droplet  diameter.  The  droplet  diameters  generated  by  the 
program  represent  the  spray  formation  characteristics  of  a  variable  area  pintle  spraybar.  The 
controlling  size  parameter  is  the  fuel  injection  pressure  drop.  These  data  are  empirical  from 
Pratt  &  Whitney  Aircraft  data.  These  five  groups  thus  represent  the  flowrate  versus  size 
distribution  unique  to  this  spraybar  type.  If  a  different  type  of  injector  is  to  be  analyzed,  the 
droplet  sizing  subroutine  must  be  rewritten. 

The  program  then  performs  a  finite  difference  solution  to  the  droplet  acceleration  and 
vaporization  equations  by  selecting  a  small  time  step  and  solving  for  the  deltas  of  liquid 
velocity,  temperature  and  a  delta  vaporized  from  the  drag  and  vaporization/heating  equations. 
The  axial  travel  of  the  liquid  droplet  during  this  time  increment  is  calculated  from  the  initial 
and  final  liquid  velocities.  The  analysis  continues  this  small  time  step  solution  until  the  axial 
distance  value  equals  the  spraybar  to  flameholder  separation  distance.  The  accelera¬ 
tion/vaporization  model  assumes  that  the  isolated  droplet  vaporization  rate  exceeds  the  rate  in 
spray  clouds  by  a  factor  of  two.  The  calculation  sums  the  amount  of  fuel  which  vaporizes  from 
all  the  droplet  size  groups  and  adds  this  to  the  amount  of  fuel  vaporized  during  injection.  The 
sum  of  these  represents  the  total  term  in  the  wake  compositional  equation. 

The  analysis  next  calculates  the  percentage  of  each  of  these  size  groups  which  is  collected 
by  impingement  onto  the  flameholder  surface.  The  collection  rate  is  caicuated  for  each  size 
group  individually  utilizing  the  droplet  diameter  which  exists  for  each  group  after  the 
vaporization  calculations.  The  total  collected  liquid  flowrate  is  summed  from  the  collection 
efficiency  for  each  group  and  the  liquid  flowrate  which  exists  in  each  group  after  the 
vaporization  calculations.  The  collecting  rate  calculation  for  each  droplet  size  group  is 
evaluated  from  droplet  trajectory  analyses  which  were  performed  and  correlated  against  the 
geometric  variables  of  the  flameholder  system.  The  details  of  this  were  presented  in  the  earlier 
Phase  I  results. 
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The  wake  recirculation  rate  is  next  evaluated,  again  following  the  earlier  reported 
analyses.  The  influence  curves  were  reduced  to  a  series  of  equations  or  curve  reading  routines 
within  the  RECIRC  subroutine.  This  subroutine  evaluates  the  aerodynamic  loading  of  the 
idealized  recirculation  zone  based  on  the  empirical  data  earlier  presented.  The  recirculation 
coefficient,  K,,  is  calculated,  and  the  loading  is  stored  for  transfer  to  the  wake  kinetics 
subroutine. 

At  this  point  the  solution  proceeds  along  the  parallel  paths  of  exercising  the  flameholder 
surface  vaporization  subroutine,  BETA  3,  and  the  wake  kinetics  subroutine. 

For  the  solution  of  BETA  3,  values  of  wake  temperature  are  assumed  for  every  200°  F 
increment  from  1000  to  5000°F.  For  each  of  these  values,  the  heat  flux  from  the  recirculation 
zone  through  the  flameholder  into  the  liquid  film  is  calculated.  The  surface  vaporization  is 
calculated  utilizing  a  20-step  finite  difference  solution  to  the  forced  convection  vaporization 
problem.  The  solution  technique  is  the  same  as  presented  earlier.  A  20-step  solution  was  found 
necessary  for  convergence  of  numerical  accuracy  and  to  assure  a  smooth  evaluation  as  the 
vaporization  rate  passes  through  to  100%.  A  check  is  made  in  the  calculation  for  the  ratio  of 
heat  flux  to  the  liquid  latent  heat  at  collection  conditions.  Whenever  the  ratio  exceeds  unity, 
the  vaporization  is  set  at  100%. 

For  each  value  of  wake  temperature,  the  calculated  vaporization  rate  is  used  in  conjunc¬ 
tion  with  the  previous  compositional  coefficients  to  define  a  wake  vapor  phase  fuel-air  ratio. 
This  temperature  versus  fuel-air  ratio  array  is  stored  for  later  use. 

The  program  transfers  next  the  WAKE  subroutine  for  the  kinetics  solution  utilizing  the 
recirculation  rate  from  RECIRC  and  the  Longwell  reactor  model  presented  earlier.  The 
subroutine  is  exercised  over  a  range  of  fuel-air  ratios.  The  fuel  air  ratios  are  started  at  0.020 
and  increased  by  0.0045.  The  solution  is  stored  in  an  array  of  efficiency  and  wake  temperature 
versus  wake  fuel-air  ratio  between  the  lean  and  rich  limits. 

These  limits  are  evaluated  at  the  level  of  aerodynamic  loading  for  the  case  being 
analyzed.  The  lean  limit  is  the  first  fuel-air  ratio  where  a  kinetic  solution  is  found.  The  rich 
limit  is  found  by  a  fine  grid  search  backwards  from  the  first  fuel-air  ratio  which  fails  to 
produce  a  solution  on  the  rich  side.  The  increased  accuracy  of  definition  of  the  rich  limit  was 
found  necessary  since  the  rich  blowout  is  the  significant  failure  mode  of  duct  flameholders. 

Once  this  array  is  generated,  the  program  searches  for  array  intersection  between  the 
BETA  3  and  WAKE  results.  If  none  is  found,  appropriate  failure  messages  are  printed  and 
execution  stops.  If  intersection  is  found,  the  convergent  results  of  the  wake  compositional 
solution  are  printed. 

The  next  subroutine,  FLAME,  performs  the  finite  difference  solution  to  the 
two-dimensional  flame  propagation  problem  presented  earlier.  The  flame  speed  base  value  is 
altered  by  the  level  of  wake  reaction  efficiency  and  by  the  percentage  of  liquid  fuel  which  is 
sloughed  from  the  flameholder  trailing  edge.  This  influence  accounts  for  two  real  effects  in  the 
augmentor  transferred  to  the  pseudo-two-dimensional  analysis.  The  wake  efficiency  is  assumed 
uniform  and  continuous  as  is  the  flame  sheet.  In  reality,  as  the  efficiency  decreases,  local 
regions  are  generated  where  flame  initiation  in  the  shear  layer  fails.  Decreasing  the  transverse 
speed  of  the  idealized  flame  sheet  is  a  method  of  describing  this  effect  to  account  for  the 
requirement  of  flame  propagation  normal  to  the  duct  into  these  unignited  regions. 
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Similarly,  the  sloughed  fuel  serves  to  quench  local  regions  of  ignition  in  the  shear  layers 
which  also  require  more  normal  flame  penetration.  The  inclusion  of  a  vaporization/sloughing 
term  accounts  for  this  in  a  two-dimensional  analysis. 

Once  the  flame  is  initiated,  the  subroutine  performs  a  finite  difference  solution  to  the 
transverse  propagation  into  the  free-stream  allowing  for  the  axial  variations  in  velocity, 
turbulence,  etc.  This  continues  until  the  exhaust  nozzle  is  reached. 

For  the  analysis  of  a  core  stream  segment  the  procedure  is  simplified  greatly  due  to  the 
complete  fuel  vaporization.  The  program  sets  the  wake  fuel-air  ratio  equal  to  the  total  fuel-air 
ratio  and  performs  the  wake  kinetics  solution  at  this  value  and  the  calculated  value  of 
recirculation.  A  gaseous  phase  turbulent  flame  is  initiated,  corrected  by  wake  efficiency,  and 
evaluated  downstream  as  before. 
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3.  Development  of  Rumble  Model  Equations 
a.  General 

The  augmentor  math  model  consists  of  a  set  of  time-dependent  equations  describing  the 
longitudinal  dynamics  of  the  flowing  air  stream  and  the  axially  distributed  combustion  process 
in  the  augmentor,  coupled  with  a  solution  technique  for  determining  stability.  These  equations 
are  linearized,  through  the  assumption  of  small  perturbations,  and  transformed  from  the 
time-domain  to  the  Laplace  transform  “S”  domain.  The  solution  technique  is  based  upon  the 
Nyquist  stability  criterion  and  consists  of  determining  whether  the  time  response  of  the 
system  to  a  small  disturbance  would  display  oscillatory  behavior  with  a  growing  amplitude. 
The  result  is  a  determination  of  stability  at  a  given  operating  point,  regions  of  operation  which 
will  cause  rumble,  and  changes  to  the  augmentor  to  make  it  rumble  free. 


b.  Development  of  Acoustic  Equations 

In  this  section  equations  are  developed  to  describe  how  velocity,  pressure,  and  density  at 
every  point  in  the  augmentor  respond  to  a  combustion  disturbance,  which  is  treated  as  a  heat 
input  to  a  flowing  inviscid  ideal  gas  stream.  Knowing  how  these  three  parameters  (velocity, 
pressure,  density)  respond  allows  calculation  of  any  other  parameter  needed,  such  as  mass 
flowrate  or  temperature.  The  first  equations  to  be  developed  are  the  three  longitudinal  wave 
equations,  which  are  applicable  between  boundaries  and  discontinuities.  Then  equations  for 
the  boundaries  and  discontinuities  are  developed.  The  wave  equations  plus  the  boundary  and 
discontinuity  equations  are  referred  to  as  the  “acoustic”  equations.  The  “combustion"  equa¬ 
tions  needed  to  complete  the  rumble  model  are  developed  in  part  c. 


Symbols  used  below  are  defined  in  the  List  of  Symbols.  For  any  section  of  augmentor 
with  rigid  walls  and  constant  cross-sectional  area,  such  as  shown  in  Figure  40,  through  which 
an  invicid  fluid  (viscosity  is  zero)  is  flowing,  the  one-dimensional  momentum,  continuity,  and 
energy  equations  are: 
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For  an  ideal  gas,  these  equations  reduce  to  the  following  non-linear  wave  equations: 
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Figure  40.  Rumble  Model  Station  Identification 


The  wave  equations  are  linearized  by  the  small  perturbation  substitutions: 

P(x,t)  =  P(x)  +  AP(x,t) 

p(x,t)  =  p(x)  +  Ap(x,t)  (89) 

C(x,t)  =  C(x)  +  AC(x,t) 

V(x,t)  =  V(x)  +  AV(x,t) 

q(x,t)  -  q(x)  +  Aq(x,t) 

Second  order  terms  are  neglected  in  making  the  substitutions. 


To  simplify  notation,  the  following  substitutions  are  made  which  normalize  the  change 
in  each  variable  by  its  steady-state  value: 


AP  AV 

P  ’  v  '  V 

The  linearized  version  of  equations  (88)  becomes: 
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Taking  the  Laplace  transform  with  respect  to  time,  with  zero  initial  conditions,  and 
letting  subscripts  1  &  2  stand  for  the  upstream  and  downstream  stations  respectively  (see 
Figure  40),  the  general  solution  to  equations  (91)  becomes: 
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In  equations  (93)  the  first  equation  describes  downstream  running  sonic  waves  of  the  form  P' 
+  V',  traveling  at  sonic  speed  plus  throughflow  velocity.  The  second  equation  describes 

upstream  running  sonic  waves  of  the  form  P'  -  7M  V',  traveling  at  sonic  speed  minus 
throughflow  velocity.  The  third  equation  describes  entropy  waves.  P'  —  yp',  drifting  down¬ 
stream  at  throughflow  velocity. 
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The  entropy  waves  become  more  apparent  from  the  expression  for  the  entropy  of  an 
ideal  gas: 

^  =  S  =  P  -  V  (94) 

The  entropy  waves  are  related  to  temperature  by : 

yr  =  S'  +  (7-I)  P  (95) 


It  is  through  equation  (95)  that  the  drifting  hot  and  cold  combustion  products,  or  entropy 
waves,  are  accounted  for  in  the  rumble  model.  Temperature  changes  produced  as  the  entropy 
waves  strike  the  exhaust  nozzle  create  waves  which  then  travel  back  upstream  at  sonic  speed. 

Equations  (9M)  are  not  useful  until  the  integrals  are  evaluated,  which  will  require 
definitions  of  V(x),  C(x),  q(x),  P(x)  and  some  assumptions  that  will  allow  integration  of  q'(x,S), 
rf|/lx,S),  d<;’(x, S),  and  (i |/(x,S).  To  complete  the  solution  the  augmentor  is  divided  into  several 
"short"  sections,  each  of  length  G,  for  each  of  which  it  can  be  assumed: 
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The  small  static  pressure  drop  in  an  augmentor  justifies  assumption  (a).  A  linear  temperature 
rise  throughout  a  section  of  length  G  is  a  good  approximation,  which  justifies  assumption  (b). 
Assumption  (c)  is  the  equation  for  a  “drifting  burning  particle”  releasing  heat  at  a  constant 
volumetric  rate  as  it  drifts  down  the  augmentor.  A  more  detailed  explanation  of  this 
assumption  will  be  provided  in  part  c  (Development  of  Combustion  Equations).  To  justify  the 
constant  steady-state  heat  release  rate  (q)  consider  the  steady-state  version  of  the  energy 
equation  (third  in  equations  (88)). 
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With  appropriate  substitutions,  the  equation  reduces  to: 


q. 

p 


/  _7_  )  R  w  dT 
'  y-1 '  P  A  dx 


V  dP 
P  dx 


Since  =  0  and  ~T~  =  constant,  then 
dx  dx 


1 

P 


=  constant  = 


7  C,M,  /  T, 
7-1  I  V  T, 


68 


(96) 


For  a  "short"  section  of  length  G  the  integration  of  f)y'(x,S)  in  equations  (93)  can  be  carried  out 
as  follows: 
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Similar  treatment  allows  integration  of  dc'(x,S)  and  liy (x,S)  in  equation  (93).  To  determine 
how  "short"  a  section  mint  he  for  the  solution  to  be  valid,  the  resulting  rumble  model  was 
exercised  repeatedly  while  decreasing  the  section  length  (by  adding  more  stations  in  the 
combustion  zone).  As  the  section  length  decreases,  the  result  will  rapidly  approach  an  exact 
solution.  It  was  found  that  section  lengths  shorter  than  about  20  inches  were  unnecessary. 


With  the  above  assumptions,  equation  (93)  becomes: 
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For  convenience  in  programming  equations  (96)  on  the  computer  the  following  identity 
substitutions  were  made: 
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where  by  definition: 
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The  time  constants  in  equations  (98)  were  evaluated  based  upon  the  steady-state 
through-flow  and  sonic  speed  profiles  created  by  the  linear  temperature  gradient. 
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C(x)  *  C, 


Then  the  time  constants  in  equations  (98)  become: 
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'This  completes  the  development  of  the  wave  equations. 


(103) 


Equations  (97)  are  applied  throughout  the  augmentor  between  any  two  stations  between 
which  there  is  no  discontinuity.  The  station  designations  used  for  the  rumble  model  are  shown 
in  Figure  40.  In  applying  the  equations,  the  general  subscripts  1  and  2  are  replaced  by  the 
actual  upstream  and  downstream  station  numbers,  respectively.  Referring  to  Figure  40,  they 
are  applied  between  stations  (1)  —  (2),  (2)  —  (3),  (4)  —  (5),  (5)  —  (10)  and  (10)  —  (11). 
Between  stations  (1)  through  (5)  and  between  stations  (10)  —  (11)  there  is  no  heat  addition, 
and  so  the  heat  addition  terms  q/P  are  set  to  zero.  The  heat  addition  terms  for  the  combustion 
zone,  stations  (5)  —  (10),  are  discussed  in  part  C. 


Discontinuities  occur  at  the  pressure  drop  locations,  stations  (2)  and  (3).  These  are 
modeled  as  small  incompressible  resistive  pressure  drops  of  zero  length.  The  continuity  and 
energy  equations  are  also  applied. 
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The  equations  are  linearized  and  normalized  as  before  to  yield: 
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In  applying  equations  (105)  to  a  given  pressure  drop  the  general  subscripts  2  and  3  are 
replaced  by  the  actual  upstream  and  downstream  station  numbers,  respectively.  For  conve¬ 
nience  in  programming,  equations  (105)  were  combined  with  the  wave  equations  (97)  to 
eliminate  the  need  for  two  stations  at  each  pressure  drop.  It  is  the  combined  equations  which 
appear  in  the  rumble  model  listing. 

A  junction  occurs  where  the  core  stream  and  fan  stream  enter  the  augmentor  and  form 
the  overall  augmentor  stream  (stations  (3),  (3H)  and  (4)  ).  Again  applying  continuity, 
momentum  and  energy: 
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For  the  Swirl  augmentor,  the  momentum  equations  at  stations  (3)  -  (4)  and  (3H)  •  (4) 
are  modified  to  account  for  the  possibility  of  different  pressure  drops  across  the  fan  and  core 
swirl  vanes.  The  linearized  version  of  the  momentum  equations  for  the  Swirl  augmentor 
becomes: 


(108) 
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Definition  of  the  upstream  and  downstream  boundary  conditions,  at  the  fan  and  at  the  noz 
zle,  respectively,  will  complete  the  acoustic  equations.  The  fan  was  assumed  to  be  delive  ring 
a  constant  mass  flowrate  through  the  fan  OD  (defined  as  that  portion  of  the  fan  between  i  he 
fan  splitter  and  fan  tip)  and  through  the  fan  ID  (defined  as  that  portion  of  the  fan  between 
the  centerline  and  the  fan  splitter).  It  was  also  assumed  that  the  temperature  of  the  fan  dis 
charge  flow  could  be  taken  as  time  invarient  (also,  because  of  the  low  Mach  number  at  fan 
discharge,  total  and  static  temperatures  can  be  used  interchangeably).  To  account  for  the 
presence  of  a  core  engine,  and  explore  any  possible  attendant  interaction  with  fan  duct 
acoustics,  a  simple  first  order  lag  representation  of  the  core  engine  was  incorporated  into  the 
rumble  model.  The  core  engine  was  represented  as  a  compressor  delivering  constant  corrected 
air  flow  (corrected  to  compressor  face  conditions)  into  a  lumped  volume.  Flow  out  of  the 
volume  exited  through  a  choked  turbine  to  emerge  at  station  (3H).  The  resulting  transfer 
function  for  the  core  engine  is: 


W,H  =  1 

Pc  1  +  T  CORE  S 

Where: 

W  ,H  =  mass  flowrate  at  station  (3H) 

P'r  =  static  pressure  at  the  compressor  face 

r cohk  =  core  engine  time  constant 


(109) 


A  default  value  of  vohk  ;  0.005  seconds  is  built  into  the  rumble  model.  A  different  .alue  can 
be  input  bv  the  user,  and  is  calculated  as  the  mass  of  air  in  the  core  engine  volume  divided  by 
the  mass  flowrate  of  air  through  the  core  engine.  Proximity  of  the  fan  splitter  to  fan  discharge 
also  affects  the  boundary  condition  at  the  fan.  Two  cases  were  considered  and  are  built  into 
the  rumble  model.  In  the  first  case,  called  the  "proximate"  splitter  configuration,  the  fan 
splitter  is  assumed  to  be  so  close  to  fan  discharge  that  no  communication  can  occur  between 
the  fan  duct  and  the  core  engine  across  the  fan  splitter.  For  this  case,  the  boundary  condition 
at  the  fan  becomes: 

Pc  W'iH  -  0 

W,  -  P  ,  +  V,-  •=  0  (110) 

T,  P,  -  p,  =  0 
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In  the  second  case,  called  the  “remote”  splitter  configuration,  the  fan  splitter  is  as¬ 
sumed  to  be  sufficiently  remote  from  fan  discharge  to  allow  perfect  communication  between 
the  fan  duct  and  the  core  engine  across  the  fan  splitter.  For  this  case,  the  boundary  condition 
at  the  fan  becomes: 
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This  completes  the  definition  of  the  upstream  boundary  condition.  It  is  of  interest  to 
note  that  entropy  waves  are  created  by  sonic  wave  reflections  at  the  upstream  boundary. 
Since  an  entropy  perturbation  is  Sj  =  P'.  -  ,  and  at  the  boundary  p’j  =  P'j ,  then 

S'  =  (1  -  y)  P'( .  A  similar  argument  will  show  that  entropy  waves  are  also  created  at  the 
pressure  drops  (stations  (2)  and  (3)  ).  These  are  automatically  accounted  for  in  the  rumble 
model,  but  are  of  minor  importance  compared  to  the  entropy  waves  created  in  the  com¬ 
bustion  zone  by  combustion  disturbances. 

The  downstream  boundary  condition  is  based  upon  the  presence  of  a  “short”  nozzle 
just  downstream  of  station  (11),  for  which: 
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P„  =  P„  /nozzle  throat  static  pressure 


When  linearized,  the  downstream  boundary  condition  becomes: 
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where: 
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It  is  also  of  interest  to  note  that  for  choked  flow, 

then  KNOZ  =  0  and: 

V,,  -"^(Pn  —  Pu)  ="2‘Tn  (114) 

substituting  from  equation  (102): 

v"  =^rSii  +1rrp"  (ll5) 

This  equation  directly  relates  how  entropy  waves,  as  well  as  pressure  disturbances,  striking  a 
choked  nozzle  will  produce  a  velocity  disturbance. 

This  completes  the  acoustic  equation  development.  These  equations  describe  the 
response  of  pressure,  velocity  and  density  throughout  the  augmentor  to  a  disturbance  in 
combustion.  Development  of  the  corresponding  combustion  equations,  which  describe  how 
combustion  throughout  the  augmentor  will  respond  to  disturbances  in  pressure,  velocity  and 
density,  is  presented  in  the  following  section. 

c.  Development  of  Combustion  Equations 

Development  of  the  combustion  equations  for  the  V-gutter  flameholder  augmentor  is 
presented  first.  Then  the  combustion  equations  for  Vorbix  and  Swirl  augmentors  are 
presented. 

For  the  V-gutter  flameholder  augmentor  two  combustion  streams,  the  fan  stream  and 
the  core  stream,  are  treated.  This  is  necessary  to  be  able  to  account  for  the  different  com¬ 
bustion  characteristics  of  the  fan  and  core  streams.  The  two  streams  can  have  different 
flameholder  designs  and  fuel-air  ratios  as  well  as  different  flameholder  approach  temperatures 
and  velocities,  causing  the  two  streams  to  have  different  efficiency  vs.  fuel-air  ratio  charac¬ 
teristics.  in  addition,  the  fan  stream  is  preceded  by  a  long  fan  duct  which  can  exhibit  lon- 
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gf.u'un  n  sonance  at  the  low  frequencies  associated  with  rumble.  The  core  stream  is  pro¬ 
ve.:  ”1  i  v  -  Tort  section  terminating  at  turbine  discharge,  which  is  much  less  responsive  at 
low  ill  qUCUClCS. 

IT.,  basic  approach  taken  for  the  rumble  model  was  to  model  combustion  disturbance;, 
in  iiu>  f.ni  and  core  streams  independently,  accounting  for  the  individual  properties  of  each 
stream.  The  resulting  two  combustion  disturbances  (calculated  as  volumetric  heat  release 
rate  disturbances)  were  then  simply  added  to  form  a  single  overall  disturbance.  The  overall 
disturbance  was  then  distributed  evenly  over  the  total  cross-sectional  area  of  the  augmentor, 
which  \<  as  taken  to  consist  of  a  single  overall  stream  with  mean  mixed  properties  This  ap¬ 
proach  accounts  for  the  different  combustion  i  haracteristics  of  the  fan  and  core  streams, 
while  u\  aiding  the  complexities  associated  with  a  rigorous  treatment  of  the  radial  as  well  as 
the  axial  distribution  of  combustion  throughout  the  augmentor. 

Kxperience  with  modeling  the  combustion  process  as  a  plane  heat  addition  with  all 
combustion  taking  place  in  zero  length  had  shown  that  the  resulting  predictions  of  rumble 
were  sensitive  to  the  axial  location  chosen  for  the  plane.  Since  combustion  actually  takes  place 
over  a  distance  of  .‘10  to  00  inches,  it  was  decided  that  the  axially  distributed  nature  of  the 
burning  should  be  accounted  for.  This  was  accomplished  hv  dividing  the  combustion  zone  into 
a  number  of  axial  sections,  each  of  length  f,  as  explained  in  part  b.  ‘‘Development  of  Acoustic 
Kqu.it  ions.” 

Combustion  equations  used  in  the  rumble  model  are  based  upon  an  extension  of 
empirical  steady-state  processes  to  the  case  of  time  variant  flow.  A  schematic  of  the 
steady  state  processes  is  shown  in  Figure  41.  Consider  first  that  the  augmentor  contains  only 
the  Ian  stream.  An  identical  set  ol  equations  will  exist  for  the  parallel  core  stream.  Following 
a  particle  of  air  as  it  moves  through  the  augmentor,  the  following  steps  will  occur: 

»  Particle  of  air  picks  up  fuel  as  it  crosses  the  spray  bar 

*  Particle  drifts  at  through  flow  velocity  to  the  flamoholuer.  station  i  1 1 

«•  Particle  is  ignited  by  the  flamehoider  wake  as  it  drifts  from  the  flame- 
holder,  to  the  beginning  of  the  combustion  zone,  station  (5)  (defined  as 
the  location  where  the  bulk  fluid  temperature  begins  to  rise  sharply). 

*•  Particle  drifts  and  burns  from  station  (5)  to  the  end  of  the  combustion 
zone,  station  (10)  (defined  as  the  location  where  bulk  fluid  temperalur- 
.  cases  its  sharp  rise ). 


!‘  "-'as  determined  (see  equation  (9(5))  that  for  a  linear  temperature  grad  it  nt.  the 
,-i caiK -state  volumetric  heat  release  rate  in  a  section  of  length  U  in  the  augmentor  could  he 
'•‘ken  ...  independent  of  axial  position.  This  implies  that  at  steady  state,  a  particle  e.l  luel-air 
mi.-.tnre  drifting  and  burning  through  a  section  of  the  augmentor,  has  a  volumetric  heat 
I'vCn-i  roe  that  is  independent  ol  axial  position.  The  rate  can  he  computed  direct  Iv  from  the 
(lu'.wab'.  ideal  temperature  rise,  efficiency,  and  combustion  zone  volume  of  the  augmentor. 
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For  tiie  fan  stream,  instantaneous  approacn  eui'ditions  are  tak  to  <>*•  <nv  'nMwut.a-.ivou» 
conditions  at  station  (3)  Because  of  tin:  large  pressure  drop  in  tne  fuel  spray  bar  injector 
changes  in  fuel  flow  m  response  to  augmontor  pressure  at  the  spray  bar  art  small  compared  to 
changes  in  air  flow.  Consequently,  fuel  flow  can  be  considered  constant,  and  the  luel-air  ratio 
of  the  particle  as  it  crosses  the  spraybar  is  determined  by  changes  in  air  flow  only. 


r ,  4  constant 

•S/H  =  W, 


A  period  of  time,  /p,  ,  is  required  f,r  tiu-  particle  to  drif'  from  th«  :.;.ra;  n-.r  «,■  the 
flainetiotiier.  Therefore,  t he  fuel-air  ratio  >1  the  particle  wticn  it  readies  the  flameholder  <-an  be 
expressed  as: 


FAc(t)  =  FAS  u  (t  -  rDC) 


(118) 


At  the  ignition  plane  (flameholder)  the  particle  has  a  “potential”  volumetric  heat  release 
rate  of : 


The  ideal  temperature  rise  is  a  function  of  the  fuel-air  ratio  of  the  particle  (effects  of  ap¬ 
proach  temperature  and  pressure  are  negligible).  The  efficiency  is  assumed  to  be  a  function  of 
the  fuel-air  ratio  and  the  approach  pressure,  temperature  and  velocity. 

Tie  -  fcn(FA,.)  (120) 

ip  =  Vfon(FAc,  P„  T„  V.) 


The  particle  crossing  the  flameholder  will  begin  burning  after  a  time  H4/V4,  which  is  the 
time  required  to  drift  from  Station  (4)  to  (5)  while  being  ignited.  When  it  begins  burning  at 
Station  (5),  the  heat  release  rate  of  the  particle  will  be  (X  =  0  at  Station  (5)  ): 

qlo.t)  =  q,  (t-#/V,) 

At  some  station,  X  distance  downstream  of  Station  (5),  the  local  heat  release  rate  will 
become  that  of  the  particle  after  an  additional  time  delay, 

t 

M-.  J  dx/V(x), 


(122) 


’I 


which  is  the  time  required  to  drift  from  Station  (5)  a  distance  X  at  through-flow  velocity 
V  (X).  Then.at  a  location  X  in  the  combustion  zone,  the  heat  release  rate  will  be: 

q(x,t)  =  q(o,t  -  rE)  023) 

The  linearized  versions  of  equations  (117)  through  (121),  written  in  terms  of  the  Laplace 
transform  of  the  normalized  variables  are: 


FAs,b  =  -  W; 
fa;  =  FA'H/8  e_Toc  s 
q;  =  w;  +  T,'  +  n' 


~  f  FA  aTt  ' 
1|c  L  T,  aFA  . 

FA; 

Jc 

•  J  FA  _aij_  1 
nc  L  v  aFA  J 

F*  +  [fS-I  Pi+[ 

solving  for  qc: 

'i  1 

q,  -  [  t-  { 

r  FA  aT,l  [  FAjaJ 

L  T,  aFAJ  t  +  L  r,  aFAJ 

.  |e‘,DCS  J 

+  [■?■?■]« 

+  [  7?r]  cTi  +  [  1T7V. 

v; 

1  c 

Xj>l 

ij  aV 


w. 


(124) 


(125) 


A  corresponding  equation  for  the  core  stream  can  be  directly  written  by  changing  sub¬ 
script  “C”  to  subscript  “H”,  and  changing  the  reference  approach  station  from  (3)  to  (3H). 


*-[  ?*].}•-•  ]w- 
’  [7*-  Lp“  +  [  +  L 


The  total  volumetric  heat  release  rate  (subscript  “T”)  is  formed  by  adding  the  heat  release 
rates  of  the  fan  and  core  streams: 

q.vT  =  Qt  =  Qc  +  Qh  =  qcvc  +  qHvH  027) 


or,  in  normalized  form : 

«-[$]«♦[$] 


Oh 


(128) 


Equation  ( 128)  computes  the  instantaneous  volumetric  heat  release  rate  of  a  particle  of 
combined  fan  stream  and  core  stream  fuel-air  mixture  when  the  particle  reaches  the 
flameholder.  The  term  “potential”  is  applied  because  the  particle  has  not  yet  been  ignited. 
The  particle  is  ignited  by  the  flameholder  wake  as  it  drifts  a  distance  P4  at  velocity  V4.  The 
particle  begins  releasing  the  “potential”  heat  at  Station  (5),  as  defined  by  equation  (122).  To 
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■  .  l,j.  :  i1  ■  ,  i,rt  has  ijiiv.  lii-tn  addc-i  '  .  ( h.  j t i.’ !, i.-i i! i r  i  ,nlv  i  lie  tun 

■..-an.  ..  •ligmallv  l*.  •..-:idrr*sl).  mi  me  t.*«ai  beat  rt  lease  <>f  both  streams  is  being 

treated  equation  1122)  is  rewritten  as: 


q(o,t)  =  qt  (t-tj  V4) 
Linearized : 


(129) 


i  •*•,:>  •  q,  (t-ijvj  (130) 

Equation  (130)  simply  adds  a  delay  into  the  system  which  allows  tailoring  the  axial  location  of 
the  beginning  of  the  combustion  zone.  For  convenience  in  programming  the  equations,  this 
delay  is  added  to  the  drift  delay  in  the  combustion  zone  (r^)  to  form  an  overall  particle  drift 
delay  from  the  flameholder. 

Ty  =  *J  V.  +  Tg  (131) 

The  particle  then  releases  heat  throughout  the  combustion  zone  as  defined  by  equation 
(123).  rhi  linearized  version  of  which  is: 

q(x,t)  =  q(o,  t-rg)  (132) 

This  equation  represents  the  augmentor  heat  release  based  on  steady-state  conditions.  To 
accurately  model  the  combustion  process  in  the  augmentor,  the  dynamics  of  the  flameholder 
wake  must  be  included.  These  dynamics  were  incorporated  in  the  heat  release  term  as  a  first 
order  log: 

_  Qoul/steadv  state  (133) 

Mout/dynamic  1  ^g 

The  dynamics  involved  in  calculating  the  flameholder  wake  time  constant  r  are  discussed  in 
Section  Il.B-2. 


Equation  (131)  was  presented  in  part  b,  (Development  of  Acoustic  Equations),  and  used 
to  evaluate  integrals  in  equation  (93).  The  combustion  equations  require  that  the  following 
information  about  the  steady-state  operating  point: 


Qc  ‘ 

'  Qh  1  r  FA 

»T,  1 

r 

FA  tr,  1 

Q,  -1 

L  Q,  -1  •  f-  Tt 

sFA  J 

C.H  *  L 

ij  *FA  J  c.h 

P  flrj 

[  lil 

1 

r  ijul  i 

t)  aP  - 

C.H  'L  V  *T 

J  C.H 

,  and 

L  i)  ftV  J  C.H 

The  heat  release  rate  ratios  Qq/Qj  and  are  computed  in  the  program  from  conditions 

known  about  each  augmentor  stream: 


Qc  _ _ (BPR  T|C  nc) _ 

Q(  (BPR  T|C  rfc)  +  (T|H  ih) 

Qh  _ _ (TiH  ith ) _ 

Qt  (BPR  T,c  ijc)  +  (T|H  >j h) 


(134) 
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are  computed  in  the  program  from  a 

C.  H 

subroutine  curvefit  of  the  ideal  temperature  rise  curve.  A  graphical  definition  of  the  term  is 

supplied  in  Figure  42.  The  partial  derivative  terms  involving  efficiency  are  computed  in  the 
F!  ''holder  Combustion  Model  and  supplied  directly  to  the  rumble  model.  Alternately,  they 
mu  e  computed  from  empirical  data  and  be  input  by  the  user.  The  graphical  definition  of 
terms  is  similar  to  that  of  Figure  42. 


The  partial  derivative  terms 
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Figure  42.  Ideal  Temperature  Rise  far  Constant  Pressure  Combustion  of 
Hydrocarbon  Fuels 


This  completes  the  combustion  equation  development  for  the  V-gutter  flameholder 
model.  All  of  the  above  equations  apply  to  the  Vorbix  and  Swirl  augmentors  except  as 
noted  below. 

For  the  Vorbix  and  Swirl  augmentors,  independent  heat  release  rates  for  the  fan  and 
core  streams  cannot  be  identified  because  of  the  flow  mixing.  In  addition,  the  effects  of 
pilot  fuel-air  ratio  on  augmentor  combustion  efficiency  must  be  accounted  for.  Equation 
(119)  is  again  applied,  but  on  an  overall  basis  only. 

Q>  -  ~  W4T,,  (135) 

The  overall  fuel-air  ratio  is  computed  from  total  mixed  air  flow  at  station  (4). 

FA  =  £2™^  (136) 

W4 
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The  overall  ideal  temperature  rise  is  a  function  of  overall  fuel-air  ratio.  The  efficiency  is  as¬ 
sumed  to  be  a  function  of  overall  fuel-air  ratio,  pilot  fuel-air  ratio  and  pressure  at  station  (4 1 


T,  =  fch(FA) 

4  =  fch(FA,FAP,P,) 


Then  for  the  Vorbix  and  Swirl  augmentors,  the  instantaneous  “potential”  volumetric  heat 
release  rate  of  a  particle  of  mixture  when  the  particle  reaches  station  (4)  is: 
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r  fa 

aT,  1 
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FA 

34 

l  T, 

aFA  J 
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aFA  JJ 

FAP 

»4 

FAR  +  [ 

34 

1  P 

>7 

aFAP 

4 

aP 

J  P* 

Equation  (139)  applies  to  both  the  Vorbix  and  Swirl  augmentors,  and  is  equivalent  to  equation 
1 1 28)  for  the  V  gutter  augmentor.  The  Vorbix  and  Swirl  augmentors  differ  in  pilot  location. 
The  Swirl  has  (he  pilot  at  fan  duet  exit,  so  that  airflow  through  the  Swirl  pilot  is  proportional 
to  fan  duct  exit  flow,  W  The  vorbix  has  the  pilot  near  midspan,  radially,  and  slightly  aft  of 
Stations  Cl)  and  Cill),  so  that  airflow  through  the  Vorbix  pilot  is  proportional  to  total  flow. 
VV,.  Then,  since  fuel  flow  into  both  pilots  is  constant: 


Swirl:  FAP  =  -  W, 
Vorbix:  FAP  =  --  W, 


(140) 


For  convenience  in  programming,  can  be  replaced  by 
W,  W,  +  W.„ 


VV  J-_MUL_1W  e 

*  L  1  *  bpr J  3 


i 


1  4  BPR 

Substituting  (109)  and  (110)  into  (108): 

FA  aT 
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(141) 
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Vorbix:  q 
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Equations  (142)  and  (142)  replace  equation  (128).  All  other  combustion  equations  are 
identical  to  those  developed  lor  the  V-gutter  flameholder  augmentor.  The  partial  derivatives  in 
equations  (142)  and  (148)  must  he  computed  from  empirical  data  and  be  input  by  the  user. 

This  completes  development  of  the  combustion  equations.  For  the  solution  technique, 
based  upon  applying  the  Nyquist  criterion  to  the  open  loop  transfer  function  (01. TF),  the 
OLTF  is  formed  by  renaming  q,’  to  q'|N  in  equation  (180)  and  by  renaming  q,  to  q'oi  I  <n 
equations  (128),  (142),  and  (142). 

d.  Phase  I  Results 

Phase  1  of  this  program  addressed  the  areas  within  the  current  model  which  had  shown 
the  need  for  improvement.  In  the  systems  dynamics  portion  of  the  model  these  areas  included: 

•  The  effect  of  liquid  fuel  droplet  dynamics  in  the  spraybar-to-flameholder 
region. 

•  The  modification  of  the  model  to  account  for  different  combustion  charac¬ 
teristics  in  the  fan  and  core  streams. 

•  The  addition  of  a  first  order  lag  on  the  heat  release  rate  due  to  flameholder 
wake  dynamics. 

These  modifications  and  their  impact  on  the  output  of  the  rumble  model  are  discussed  in 
this  section. 


(1)  Liquid  Fuel  Droplet  Dynamics  in  Spraybar-to-Flameholder  Region 

The  rumble  model,  as  developed  ur  the  I.o-Frequency  Augmentor  Instability  In¬ 
vestigation.  predicted  the  first  mode  of  rui-n.ie  at  higher  frequencies  than  experienced  during 
the  FSEK  tests.  Also,  a  second  or  third  mode  of  rumble  was  predicted  that  was  not  observed 
during  the  FSF.R  tests  (Figure  42).  One  possible  explanation  is  that  there  is  a  lag  in  the  heat 
release  rate  (Q  out)  that  is  not  accounted  for  in  the  model.  The  model  assumes  gaseous  fuel 
particles  that  move  in  phase  with  the  air  How  between  (he  spraybars  and  the  flameholders.  In 
fact,  liquid  particles  exist  in  this  region  and  should  exhibit  a  lag  in  directional  change  which 
would  cause  a  lag  in  the  heat  release  rate.  During  the  I.o-Frequency  Augmentor  Instability 
Investigation  a  lug  was  incorporated  on  the  heat  release  rate  (Q  out)  to  verify  the  assumption. 
The  rumble  model  was  exercised  at  several  altitude  test  points  with  lags  ranging  from  0.001  to 
0.008  sec.  A  lag  of  0.00.8  sec  indicated  best  agreement  with  test  results.  Figure  44  presents  the 
same  test  point  as  presented  in  Figure  48  with  a  0.002  sec  lag  on  the  heat  release  rate.  The 
second  and  third  rumble  modes  are  not  predicted  and  the  first  mode  of  rumble  at  52  Hz  is  in 
perfect  agreement  with  FSF.R  test  results. 
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During  Phase  I  of  this  contract,  an  analytical  investigation  was  conducted  into  the 
physics  of  the  liquid  fuel  particle  dynamics  in  the  spray bar/flameholder  region.  The  augmentor 
fuel  system  model  was  used  to  predict  liquid  fuel  particle  sizes  for  various  augmentor/flight 
conditions.  The  droplet  dynamics  were  modeled  by  applying  the  equation  F  =  ma,  where  F  is 
the  drag  force  on  the  droplet,  m  is  the  mass  of  the  droplet  and  a  is  the  acceleration  of  the 
droplet.  This  results  in  a  first  order  differential  equation  for  the  velocity  of  the  fuel  droplet: 


dVf 

dt 


+ 


Td 


(144) 


where  Vf  is  the  fuel  droplet  velocity,  Va  is  the  velocity  of  the  airstream,  and  rd  is  the  fuel 
droplet  time  constant  which  is  a  function  of  the  droplet  diameter  and  the  drag  coefficient, 
rd  =  f  (d,  Cp).  In  the  frequency  domain,  using  normalized  variables 


I  ) 


A  (  ) 
(  I 


(145) 


the  fuel  droplet  velocity  is  given  by: 


'  1  +  rdS  •  (146) 

The  fuel  air  ratio  is  given  by: 

PA  -  (UV 

where  Wf  and  W.,  are  the  fuel  flow  rate  and  air  flow  rate,  respectively.  Linearized,  this  becomes 
FA'  =  W(  -  Wa'  where 

W'f  =  Of  +  Vf  +  0'  and  Wa'  =  pa’  +  Va-  +  '0  (148) 

and  pi  and  pa  are  the  fuel  and  air  densities. 

Making  these  substitutions,  the  equation  for  the  fuel-air  ratio  becomes: 

FA'  =  TT^i"  "  "  v°' or  FA'  =  7 T^s  Va'  '  Pa  (149) 

This  is  the  equation  that  was  added  to  the  rumble  model. 

In  the  original  model,  fuel  flow  was  assumed  constant  and  the  change  in  fuel-air  was 
given  by  FA'  =  -  Va.  -  pa . 

The  effect  of  droplet  dynamics  is  to  modify  the  effect  in  air  velocity  by  the  term 


TdS 

1  +  T<jS 

Physically,  this  term  is  a  high  pass  filter  and  attenuates  frequencies  below 


(150) 


1 

2irrH 
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(151) 


Results  of  the  rumble  model  with  the  addition  of  fuel  droplet  dynamics  (r^  =  0.0053  sec) 
are  shown  in  Figure  45.  Instabilities  occur  at  75  Hz,  147  Hz  and  232  Hz,  which  do  not  agree 
with  engine  data.  The  combustion  lag  shown  in  the  figure  acts  as  a  low  pass  filter.  Therefore, 
the  addition  of  fuel  droplet  dynamics,  which  acts  as  a  high  pass  filter,  was  not  expected  to  give 
results  in  agreement  with  test  data.  The  effect  of  adding  both  fuel  droplet  dynamics  and  the 
0.003  sec  combustion  lag  was  also  evaluated.  Figure  46  shows  that  this  combination  is  stable. 

These  results  indicate  that  the  model,  with  the  addition  of  fuel  droplet  dynamics,  does 
not  agree  with  engine  test  data.  Therefore,  fuel  droplet  dynamics  can  be  eliminated  as  a 
critical  factor  in  the  prediction  of  rumble. 

(2)  Account  for  Different  Combustion  Characteristics  in  the  Fan  and  Core  Streams 

When  the  rumble  model  was  extended  to  account  for  the  fan  and  core  stream  of  a 
turbofan  engine,  certain  assumptions  were  made  to  simplify  the  wave  equations.  One  of  these 
assumptions  was  that  downstream  of  the  flameholder,  the  fan  and  core  streams  are  mixed  and 
come  to  a  common  temperature  at  Station  4  (Figure  47).  The  temperature  gradient  was 
assumed  to  be  linear  throughout  the  augmentor.  The  overall  heat  release  of  the  combined  fan 
and  core  streams  was  assumed  to  begin  at  Station  5  and  end  at  Station  10. 

During  Phase  I  of  this  contract,  the  rumble  model  was  extended  to  account  for  the  fact 
that  the  fan  stream  and  core  stream  have  different  axial  temperature  profiles.  Figure  48  shows 
typical  augmentor  temperature  profiles  for  various  augmentor  fuel/air  ratios.  Since  it  is 
required  in  the  rumble  model  that  the  temperature  rise  is  linear  between  two  sections  in  the 
augmentor,  the  actual  temperature  profile  had  to  be  modified.  The  temperature  profile  was 
broken  up  into  a  series  of  linear  gradients  which  approximated  the  actual  temperature  profile. 
It  was  found  that  a  profile  formed  by  combining  a  core  stream  temperature  rise  occurring  over 
the  first  three  sections  of  the  augmentor  and  a  fan  stream  temperature  rise  occurring  over  the 
last  four  sections  gave  the  best  approximation.  Figure  49  shows  the  individual  fan  and  core 
stream  temperature  profiles  used  in  the  rumble  model,  and  Figure  50  shows  the  combined 
profile.  This  improved  temperature  profile  reflects  the  different  combustion  characteristics  of 
the  fan  and  core  streams. 


Data  from  FX217-18  engine  tests  at  AEDC  were  used  to  verify  model  changes.  The  flight 
point  and  flameholder  design  selected  were  0.8/30K  and  Frithorn  flameholder  (Run  70/1100). 
Engine  data  indicated  that  this  was  a  stable  operating  point. 

The  original  rumble  model  (see  AFAPL-TR-78-82)  and  the  rumble  model  with  the 
imnroved  augmentor  temperature  profile  were  both  run  using  FX217-18  data  from  Run 
70/1 100  at  0.8/39K.  The  combustion  model  was  executed  using  a  single  streamtube  analysis  to 
generate  input  for  the  rumble  model.  Figure  51  gives  the  output,  of  the  original  rumble  model. 
An  instability  is  predicted  at  57  Hz.  Engine  test  data  indicates  stable  operation  at  this  point 
so  the  prediction  of  the  original  rumble  model  is  incorrect. 

The  output  of  the  rumble  model  with  the  improved  augmentor  temperature  profile  is 
shown  in  Figure  52. 

The  plot  indicates  an  instability  at  60  Hz.  This  is  not  in  agreement  with  test  data. 
However,  since  this  improved  temperature  profile  is  a  better  model  of  the  actual  combustion 
process,  it  was  incorporated  into  the  model  as  a  permanent  change. 
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Figure  47.  Original  Model  of  Augmcntor  Temperature  Rise 
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Figure  51.  Original  Rumble  Model.  Frithorn  Flameholder,  Single  Stream- 
tube  Analysis,  0.8/39K 


0.0  25.0  50.0  75.0  100.0  125.0  150.0  175.0  200.0  225.0  250.0 

Frequency  -  Hertz 

FD  230871 

Figure  52.  Rumble  Model  With  Improved  Augmentor  Temperature  Profile. 

Frithorn  Flameholder,  Single  Streamtube  Analysis,  0.8/39K 


(3)  Lag  on  Heat  Release  Rate  (Qoui) 

Since  the  addition  of  fuel  droplet  dynamics  in  the  spraybar-to  flameholder  region  did  not 
result  in  the  appropriate  lag  on  the  heat  release  rate,  another  source  of  this  lag  was  sought.  It 
was  found  that  there  is  a  lag  associated  with  the  dynamics  of  the  flameholder  wake  region.  The 
dynamics  involved  in  the  calculation  of  this  flameholder  wake  time  constant  are  discussed  in 
detail  in  Section  Il.B-2.  This  lag  was  incorporated  into  the  rumble  model  as  a  single  first  order 
lag  on  the  total  (duct  plus  core)  heat  release  rate  based  on  steady-state  conditions: 


Q..UI  tknuntu- 


Q,„ 


steady  stair 


+  7-S 


Where  r  is  the  flameholder  wake  time  constant. 


(152) 


The  individual  time  constants  for  each  streamtube  in  the  duct  and  core  s'  ’•earns  were 
averaged  to  give  a  single  overall  time  constant.  Typical  values  of  this  time  constant  are  0,00.) 
to  0.005  sec. 

Figure  5,'i  shows  the  output  of  the  rumble  model  for  the  Frithorn  flameholder  at  <).8/.'i9K 
with  the  addition  of  the  improved  augmentor  temperature  gradient  and  the  lag  on  the  heat 
release  rate.  Again,  the  combustion  model  was  run  using  a  single  streamtube  analysis.  The 
flameholder  wake  time  constant  calculated  for  this  case  is  0.005  sec.  The  frequency  has  shifted 
to  lower  values  and  the  gain  at  the  high  frequencies  has  been  reduced.  These  arc  the  expected 
results  of  the  addition  of  a  first  order  lag  on  the  heat  release  rate.  The  output  indicates  that 
this  is  a  stable  condition  which  agrees  with  engine  test  data. 

e.  Refined  Mode I  Evaluation 


Phase  II  of  the  current  program  concentrated  on  comparison  of  the  refined  rumble  model 
with  available  engine  test  data  taken  on  the  FlOOCi).  These  data  were  measured  at  the  AKDC 
test  (ells  on  engine  FX217-18  with  two  different  flameholder  configurations,  identified  as 
Flood)  Bill  of  Material  and  a  revised  design  called  "Frithorn.”  The  Phase  I  results  presented 
the  refined  rumble  model  output  at  a  flight  point  of  0.8  Mach  number  and  29.000  ft  altitude 
I0.8/39K),  using  a  single  fan  and  single  core  streamtube.  The  Phase  II  results  utilize  multiple 
streamtubes  to  evaluate  fuel-air  ratio  distribution.  Combinations  of  0.8/39K  and  0.8/45K  flight 
points  and  the  two  flameholder  configurations  wrere  evaluated. 

(1)  Muttistreamtube  Analysis 

The  combustion  model  was  executed  using  a  multistreamtube  analysis  in  the  core  and 
duel  for  the  Frithorn  flameholder  at  the  flight  point  0.8/39K.  This  analysis  allowed  the  user  to 
model  the  radial  and  circumferential  variations  in  geometry  and  aerod-.  'amics.  The  analysis 
consisted  of  48  fan  duct  streamtubes  and  84  core  streamtubes.  The  results  of  the  improved 
rumble  model  using  the  multistreamtube  combustion  model  output  are  shown  in  Figure  54. 
The  model  indicates  that  this  is  a  stable  operating  point  which  agrees  with  test  data.  The  gain 
in  the  heat  release  rate  is  slightly  lower  than  the  gain  calculated  using  a  single  streamtube 
combustion  model.  Since  the  fuel-air  ratios  vary  significantly  in  any  engine,  it  is  recommended 
that  the  multistreamtube  analysis  be  used. 
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Figure  53  Rumble  Model  With  Improved  Augmentor  Temperature  Gra¬ 
dient  and  Lag  on  Qoui  —  Frithorn  Flameholder,  0.8/39K,  Single 
Streamtube  Analysis 
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Iif’urc  5-1.  Humble  Model  With  Improved  Augmentor  Temperature  (!ra- 
dient  and  Lag  on  Q„ul  —  Frithorn  Flameholder,  0.8/39K ,  Multi- 
streamtube  Analysis 


(2)  Analysis  of  Additional  Data 

For  further  verification  the  model  was  executed  using  data  from  various  other  combina¬ 
tions  of  flameholder/flight  conditions. 

Using  a  multistreamtube  combustion  analysis,  the  model  predicted  rumble  at  48  Hz  for 
the  Frithorn  flameholder  at  0.8/45K  (see  Figure  55).  FX217-18  engine  data  indicated  a  rumhle 
blowout  at  0.8/4 IK  (Figure  56).  The  spectrum  analysis  at  the  bottom  of  Figure  56  indicates 
that  the  rumble  frequency  is  approximately  50  Hz. 


The  gain  calculated  in  'he  rumble  model  at  0.8/39K  is  0.6  and  the  gain  at  0.8/45K  is  1.4. 
Somewhere  between  these  two  points  the  model  would  predict  the  transition  from  stable 
operation  to  rumble.  Assuming  a  linear  variation  in  gain  from  0.8/39K  to  0.8/45K,  we  expect 
the  model  to  predict  rumble  starting  at  approximately  0.8/42K.  The  altitude  predicted  for  the 
onset  of  rumble  for  the  Frithorn  flameholder  is  in  good  agreement  with  test  data.  The  rumble 
frequency  predicted  by  the  model  also  compares  well  with  engine  data.  Table  I  gives  a 
summary  of  the  results  for  the  Frithorn  flameholder. 

The  model  was  run  for  the  B/M  flameholder  at  0.8/39K.  Figure  57  shows  the  rumble 
model  output  for  this  case.  The  model  predicts  rumble  (gain  ~-2)  at  48  Hz.  When  execution  of 
the  rumble  model  was  attempted  at  0.8/45K.  an  error  was  encountered.  The  combustion  model 
indicated  that  the  rich  limit  on  fuel -air  ratio  in  the  duct  was  exceeded.  Therefore,  the  rumble 
model  could  not  be  run  at  this  flight  condition. 

Engine  data  from  FX217-18  for  the  B/M  flameholder  at  0.8/39K  indicates  stable 
operation.  A  rumhle  blowout  occurred  at  0.8/49K  (see  Figure  58).  The  spectrum  analysis  shows 
the  frequency  of  rumble  to  be  approximately  50  Hz. 

Rumble  model  predictions  for  the  B/M  flameholder  do  not  agree  with  FX217-18  test 
data.  The  model  predicts  rumble  at  0.8/39K  while  engine  data  indicates  stable  operation.  The 
combustion  model  is  predicting  a  value  of  -7  for  ZEFC 

(■T'SFA  >duct) 

for  this  case.  ZEFC  is  a  main  factor  involved  in  the  prediction  of  rumble.  This  large  negative 
value  of  ZEFC  is  the  reason  for  the  model  prediction  rumble  at  0.8/39K  for  the  B/M 
flameholder.  Also,  the  combustion  model  would  not  generate  any  output  for  the  flight  point 
0.8/45K.  The  engine  was  operating  at  this  point.  Therefore,  the  combustion  model  results  do 
not  agree  with  engine  data  for  the  B/M  flameholder. 
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Figure  55.  Rumble  Model  —  Frithorn  Flameholder ,  0.8/45K,  Multistreom- 
tube  Analysis 
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TABLE  I.  FRITHORN  FLAMEHOLDER 

Altitude  for  Rumble 
at  0.8  Mach  No.  Rumble  Frequency 

_ (JO _ (Hz) 

Engine  Data  —  FX217-18  41K  50 
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Figure  57.  Rumble  Model  —  B/M  Flameholder,  0.8/39K,  Mult istre amt ube 
Analysis 
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58.  FX217-18,  Run  69/1090.  R/M  Flameholder,  0.8/48K  Augmentor  Mistier  AK03 


C.  CONCLUSIONS 


The  Augmentor  Stability  Management  Program  has  resulted  in  an  extended  and  im¬ 
proved  design  system  for  turbofan  augmentors.  The  program  is  capable  of  design  evaluation 
for  performance  or  stability.  The  results  of  the  various  submodels  are  in  excellent  agreement 
with  experimental  data  from  several  sources. 

The  combustion  model  uses  a  unique  solution  for  stability  and  efficiency  which  explains 
the  relative  effects  of  fuel-air  ratio,  geometry  and  inlet  conditions.  The  model  has  shown 
reasonable  correlation  with  full  scale  engine  test  data.  The  full  program  is  a  viable  design 
system  for  engine  design  and  development. 

Perhaps  the  most  important  conclusions  are  drawn  from  the  sensitivity  of  the  results  to 
variations  in  inlet  conditions.  To  design  a  stable,  high  efficiency  turbofan  augmentor  fuel 
system  and  flameholder,  the  designer  must  have  accurate  quantitative  definition  of  the  inlet 
profiles.  This  information  should  be  considered  an  early  priority  item  in  an  engine  develop¬ 
ment  program. 
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